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Abstract 

We discuss the asymptotic behaviour of models of lattice polygons, mainly on 
the square lattice. In particular, we focus on limiting area laws in the uniform 
perimeter ensemble where, for fixed perimeter, each polygon of a given area occurs 
with the same probability. We relate limit distributions to the scaling behaviour of the 
associated perimeter and area generating functions, thereby providing a geometric 
interpretation of scaling functions. To a major extent, this article is a pedagogic 
review of known results. 



1 Introduction 

For a given combinatorial class of objects, such as polygons or polyhedra, the most basic 
question concerns the number of objects of a given size (always assumed to be finite), or 
an asymptotic estimate thereof. Informally stated, in this overview we will analyse the 
refined question: 

What does a typical object look like? 



In contrast to the combinatorial question about the number of objects of a given size, the 
latter question is of a probabilistic nature. For counting parameters in addition to object 
size, one asks for their (asymptotic) probability law. To give this question a meaning, 
an underlying ensemble has to be specified. The simplest choice is the uniform ensemble, 
where each object of a given size occurs with equal probability. 

For self-avoiding polygons on the square lattice, size may be the number of edges of the 
polygon, and an additional counting parameter may be the area enclosed by the polygon. 
We will call this ensemble the fixed perimeter ensemble. For the uniform fixed perimeter 
ensemble, one assumes that, for a fixed number of edges, each polygon occurs with the 
same probability. Another ensemble, which we will call the fixed area ensemble, is obtained 
with size being the polygon area, and the number of edges being an additional counting 
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parameter. For the uniform fixed area ensemble, one assumes that, for fixed area, each 
polygon occurs with the same probability. 

To be specific, let p m) n denote the number of square lattice self-avoiding polygons of 
half-perimeter m and area n. Discrete random variables X m of area in the uniform fixed 
perimeter ensemble and of perimeter Y n in the uniform fixed area ensemble are defined by 



F(X m = n)= , P(y n = m 



n Pm,n 2—imP m > n 

We are interested in an asymptotic description of these probability laws, in the limit of 
infinite object size. 

In statistical physics, certain non-uniform ensembles are important. For fixed object 
size, the probability of an object with value n of the counting parameter (such as the area 
of a polygon) may be proportional to a n , for some non- negative parameter a = e~@ E of 
non-uniformity. Here E is the energy of the object, and (3 = l/(fceT), where T is the 
temperature, and ks denotes Boltzmann's constant. A qualitative change in the behaviour 
of typical objects may then be reflected in a qualitative change in the probability law of 
the counting parameter w.r.t. a. Such a change is an indication of a phase transition, i.e., 
a non-analyticity in the free energy of the corresponding ensemble. 

For self-avoiding polygons in the fixed perimeter ensemble, let q denote the parameter 
of non-uniformity, 

nX m (q) =n) = J? m f 

Polygons of large area are suppressed in probability for small values of q, such that one 
expects a typical self-avoiding polygon to closely resemble a branched polymer. Likewise, 
for large values of q, a typical polygon is expected to be inflated, closely resembling a ball 
(or square) shape. Let us define the ball-shaped phase by the condition that the mean 
area of a polygon grows quadratically with its perimeter. The ball-shaped phase occurs 
for q > 1 [31]. Linear growth of the mean area w.r.t. perimeter is expected to occur 
for all values < q < 1. This phase called the branched polymer phase. Of particular 
interest is the point q = 1, at which a phase transition occurs [31]. This transition is called 
a collapse transition. Similar considerations apply for self-avoiding polygons in the fixed 
ensemble, 

Fm,n Jj 



¥(Y n (x) = m) 



Em Pm,nX r 

with parameter of non-uniformity x, where < x < oo. 

For a given model, these effects may be studied using data from exact or Monte-Carlo 
enumeration and series extrapolation techniques. Sometimes, the underlying model is 
exactly solvable, i.e., it obeys a combinatorial decomposition, which leads to a recursion 
for the counting parameter. In that case, its (asymptotic) behaviour may be extracted 
from the recurrence. 

A convenient tool is generating functions. The combinatorial information about the 
number of objects of a given size is coded in a one-variable (ordinary) generating function, 
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typically of positive and finite radius of convergence. Given the generating function of the 
counting problem, the asymptotic behaviour of its coefficients can be inferred from the 
leading singular behaviour of the generating function. This is determined by the location 
and nature of the singularity of the generating function closest to the origin. There are 
elaborate techniques for studying this behaviour exactly [37J or numerically [43J. 

The case of additional counting parameters leads to a multivariate generating function. 
For self-avoiding polygons, the half-perimeter and area generating function is 



m n 



For a fixed value of a non- uniformity parameter qo, where < go < 1, let xq be the radius 
of convergence of P(x, qo)- The asymptotic law of the counting parameter is encoded in the 
singular behaviour of the generating function P(x, q) about (xq, qo). If locally about (xq, qo) 
the nature of the singularity of P(x, q) does not change, then distributions are expected 
to be concentrated, with a Gaussian limit law. This corresponds to the physical intuition 
that fluctuations of macroscopic quantities are asymptotically negligible away from phase 
transition points. If the nature of the singularity does change locally, we expect non- 
concentrated distributions, resulting in non-Gaussian limit laws. This is expected to be 
the case at phase transition points. 

Qualitative information about the singularity structure is given by the singularity di- 
agram (also called the phase diagram). It displays the region of convergence of the two- 
variable generating function, i.e., the set of points (x, q) in the closed upper right quadrant 
of the plane, such that the generating function P(x, q) converges. The set of boundary 
points with positive coordinates is a set of singular points of P(x,q), called the critical 
curve. See Figure [1] for a sketch of the singularity diagram of a typical polygon model 
such as self-avoiding polygons, counted by half-perimeter and area, with generating func- 
tion P(x, q) as above. There appear two lines of singularities, which intersect at the point 



Xc 



Figure 1: Singularity diagram of a typical polygon model counted by half-perimeter and area, 
with x conjugate to half-perimeter and q conjugate to area. 

(x, q) = (x c , 1). Here x c is the radius of convergence of the half-perimeter generating func- 
tion P(x, 1), also called the critical point. The nature of a singularity does not change 
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along each of the two lines, and the intersection point (x, q) = (x c , 1) of the two lines is a 
phase transition point. For < q < 1 fixed, denote by x c (q) the radius of convergence of 
P(x, q). The branched polymer phase for the fixed perimeter ensemble < q < 1 (and also 
for the corresponding fixed area ensemble) is asymptotically described by the singularity 
of P(x, q) about (x c (q), q). In the ball-shaped phase q > 1 of the fixed perimeter ensemble, 
the (ordinary) generating function does not seem the right object to study, since it has zero 
radius of convergence for fixed q > 1. The singularity of P(x, q) about (x, 1) describes, for 
< x < x c , a ball-shaped phase in the fixed area ensemble, with a finite average size of a 



For points (x, q) within the region of convergence, both x and y positive, the generating 
function P(x, q) is finite and positive. Thus, such points may be interpreted as parameters 
in a mixed infinite ensemble 



The limiting law of the counting parameter in the fixed area or fixed perimeter ensem- 
ble can be extracted from the leading singular behaviour of the two-variable generating 
function. There are two different approaches to the problem. The first one consists in 
analysing, for fixed non- uniformity parameter a, the singular behaviour of the remaining 
one-parameter generating function and its derivatives w.r.t. a. This method is also called 
the method of moments. It can be successfully applied in the fixed perimeter ensemble at 
the phase transition point. Typically, this results in non-concentrated distributions. 

The second approach derives an asymptotic approximation of the two- variable generat- 
ing function. Away from a phase transition point, such an approximation can be obtained 
for some classes of models, typically resulting in concentrated distributions, with a Gaus- 
sian law for the centred and normalised random variable. However, it is usually difficult to 
extract such information at a phase transition point. The theory of tricritical scaling seeks 
to fill this gap, by suggesting and justifying a particular ansatz for an approximation using 
scaling functions. Knowledge of the approximation may imply knowledge of the quantities 
analysed in the first approach. 

In the following, we give an overview of these two approaches. For the first approach, 
summarised by the title limit distributions, there are a number of rigorous results, which 
we will discuss. The second approach, summarised by the title scaling functions, is less 
developed. For that reason, our presentation will be more descriptive, stating important 
open questions. We will stress connections between the two approaches, thereby providing 
a probabilistic interpretation of scaling functions in terms of limit distributions. 

2 Polygon models and generating functions 

Models of polygons, polyominoes or polyhedra have been studied intensively on the square 
and cubic lattices. It is believed that the leading asymptotic behaviour of such models, 
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such as the type of limit distribution or critical exponents, is independent of the underlying 
lattice. 

In two dimensions, a number of models of square lattice polygons have been enumerated 
according to perimeter and area and other parameters, see [7J for a review of models with an 
exact solution. The majority of such models has an algebraic perimeter generating function. 
We mention prudent polygons [961 [22l [8] as a notable exception. Of particular importance 
for polygon models is the fixed perimeter ensemble, since it models two-dimensional vesicle 
collapse. Another important ensemble is the fixed area ensemble, which serves as a model of 
ring polymers. The fixed area ensemble may also describe percolation and cluster growth. 
For example, staircase polygons are models of directed compact percolation [26| [28j [291 
1271 [T2| [57] . This may be compared to the exactly solvable case of percolation on a tree 
|42j . The model of self-avoiding polygons is conjectured to describe the hull of critical 
percolation clusters [60] . 

In addition to perimeter, other counting parameters have been studied, such as width 
and height, generalisations of area [SU], radius of gyration [521 El], number of nearest- 
neighbour interactions [1], last column height [7J, and site perimeter [201 EH]- Also, mo- 
tivated by applications in chemistry, symmetry subclasses of polygon models have been 
analysed [631 E2l HQl ES]- Whereas this gives rise to a number of different ensembles, only 
a few of them have been asymptotically studied. Not all of them display phase transitions. 

In three dimensions, models of polyhedra on the cubic lattice have been enumerated 
according to perimeter, surface area and volume, see [711 H021 E] and the discussion in 
section 13.91 Various ensembles may be defined, such as the fixed surface area ensemble 
and the fixed volume ensemble. The fixed surface area ensemble serves as a model of 
three-dimensional vesicle collapse [104J. 

In this chapter, we will consider models of square lattice polygons, counted by half- 
perimeter and area. Let p m ,n denote the (finite) number of such polygons of half-perimeter 
m and area n. The numbers p m .n will always satisfy the following assumption. 

Assumption 1. Form,n e No, let non-negative integers p m . n G No be given. The numbers 
Pm,n are assumed to satisfy the following properties. 

i) There exist positive constants A, B > such thatp m ^ n = if n < Am or if n > Bm 2 . 

ii) The sequence {^2 n Pm,n) m&no has infinitely many positive elements and grows at most 
exponentially. 

Remarks, i) A sequence (a n ) ng N is said to grow at most exponentially, if there are positive 
constants C, fi such that \a n \ < Cfi n for all n. 

ii) Condition i) reflects the geometric constraint that the area of a polygon grows at most 
quadratically and at least linearly with its perimeter. For self-avoiding polygons, we have 
n > m — 1. Since p m ,n = if m < 2, we may choose A = 1/3. Since n < m 2 /4 for 
self-avoiding polygons, we may choose 5 = 1/3. Condition ii) is a natural condition on 
the growth of the number of polygons of a given perimeter. For self-avoiding polygons, we 
may choose C = 1 and /i = 16. 
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Hi) For models with counting parameters different from area, or for models in higher 
dimensions, a modified assumption holds, with the growth condition i) being replaced by 
n < Am k ° and n > Bm kl , for appropriate values of ko and k%. Counting parameters 
statisfying p m , n = for n > Bm k are called rank k parameters |25j . 

The above assumption imposes restrictions on the generating function of the numbers 
p m , n . These explain the qualitative form of the singularity diagram Figure HJ 

Proposition 1. For numbers p m ,n, let Assumption^ be satisfied. Then, the generating 
function P(x, q) = J2mnP-m,nX m q n has the following properties. 

i) The generating function P(x, q) satisfies for k G N 

k / o \ k / o \ 2k 

| P(x,q)^B k (x—) P(x,q), 







dx 



^ ( x— ] P(x,q)< [q— 







dq 



where denotes coefficient-wise domination. 

ii) The evaluation P(x, 1) is a power series with radius of convergence x c , where < 
x c < 1. 

Hi) The generating function P(x,q) diverges, if x ^ and \q\ > 1. It converges, if\q\ < 1 
and \x\ < x c q~ A . In particular, for fc e No, the evaluations 



Qk 

dx k 



P{x,q) 



are power series with radius of convergence 1. 
iv) For k G N 0; the evaluations 



d k 
dq 1 



■P(x,q) 



q=l 



are power series with radius of convergence x c . They satisfy, for \x\ < x c , 

Qk 



dq' 



: P(x,q) 



Qk 

= lim -^-rP(x,q). 
, aq K 
i- 1 -1<?<1 ^ 



sketch. The domination formula follows immediately from condition i). The existence of 
the evaluations at q — 1 and formal power series also follows from condition i). 

Condition ii) ensures that < x c < 1 for the radius of convergence of P(x, 1). Equality 
of the radii of convergence for the derivatives follows from condition i) by elementary 
estimates. The claimed analytic properties of P(x,q) follow from conditions i) and ii) by 
elementary estimates. The claimed left-continuity of the derivatives in iv) is implied by 
Abel's continuity theorem for real power series. □ □ 
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Remarks, i) Proposition [T] implies that the critical curve x c (q) satisfies for < q < 1 the 
estimate x c (q) > x c q~ A . For self-avoiding polygons, the critical curve x c (q) is continuous 
for < q < 1. This follows from a certain supermultiplicative inequality for the numbers 
Pm,n by convexity arguments [38] . 

ii) Of central importance in the sequel will be the power series 

(1) 

9=1 

They are called factorial moment generating functions, for reasons which will become clear 
later. 

We continue studying analytic properties of the factorial moment generating functions. 
In the following, the notation i/i denotes the limit x — > x for sequences (x n ) satisfying 
\x n \ < x . The notation f(x) ~ g(x) as x S x means that g (x) ^ in a left neigbourhood 
of Xq and that lim x/ » xo f(x)/g(x) = 1. Likewise, a m ~ b m as m — > oo for sequences 
(a m ),(6 m ) means that 6 m 7^ for almost all m and lim ni _ j . 0O a m /b m = 1. The following 
lemma is a standard result. 

Lemma 1. Le£ (a m )meN be a sequence of real numbers, which asymptotically satisfy 

a m ~ ArT™™ 7 " 1 (m -»• 00), (2) 

for real numbers A, x c , 7, where A 7^ and x c > 0. 

T/ien, i/ie generating function g(x) = Y^=o a mX m has radius of convergence x c . If 
7 ^ {0,-1, —2, . . .}, then there exists a power series g^ re9 \x) with radius of convergence 
strictly larger than x c , such that g(x) satisfies 

{g(x) - g<™\x)) ~ ,/ r( /\ 7 (x / x e ), (3) 
(1 — x/x c y 

where T(z) denotes the Gamma function. □ 

Remarks, i) The above lemma can be proved using the analytic properties of the polylog 
function |32j. If 7 e {0, —1, —2, . . .}, an asymptotic form similar to Eq. ([3]) is valid, which 
involves logarithms. 

ii) The function the above lemma is not unique. For example, if 7 > 0, any 

polynomial in x may be chosen. We demand = in that case. If 7 < and 

g( reg \x) is restricted to be a polynomial, it is uniquely defined. If — 1 < 7 < 0, we 
have g( re9 \x) = g(x c ). In the general case, the polynomial has degree [ — 7J, compare 
[3^] . In the following, we will demand uniqueness by the above choice. The power series 
g( sm 9)(x) := (g(x) — g( re9 \x)) is then called the singular part of g(x). 

Conversely, let a power series g(x) with radius of convergence x c be given. In order to 
conclude from Eq. ([3D the behaviour Eq. (|2J), certain additional analyticity assumptions 
on g(x) have to be satisfied. To this end, a function g(x) is called A(x c ,r],(f>) -regular (or 



7 



simply A-regular) [30], if there is a positive real number x c > 0, such that g(x) is analytic 
in the indented disc A(x c ,rj,4>) := {z G C : \z\ < x c + r], \Aig(z — x c )\ > <p}, for some 
i] > and some 0, where < <f> < n/2. Note that x c ^ A, where we adopt the convention 
Arg(0) = 0. The point x = x c is the only point for \x\ < x c , where g(x) may possess a 
singularity. 



~ — 7 (x — >■ x c m A). 



Lemma 2 Q35J). Lei the function g(x) be A-regular and assume that 

1 

(1 - x/x c )~ 
if 7 ^ {0, —1, —2, . . we then have 

[x m ]g(x) ~ ^-yx^m 7 - 1 (m -> oo), 

where [x m ]g(x) denotes the Taylor coefficient of g(x) of order m about x = 0. □ 

Remarks, i) Note that the coefficients of the function f(x) = (1 — x/x c )~ J with real 
exponent 7 ^ {0, —1, —2, . . .} satisfy 

[x m ]f(x) ~ f^y^c^ 7 " 1 ( m -» oo)- ( 4 ) 

This may be seen by an application of the binomial series and Stirling's formula. For 
functions g(x) ~ /(^), the assumption of A-regularity for g(x) ensures that the same 
asymptotic estimate holds for the coefficients of g(x). 

ii) Theorems of the above type are called transfer theorems [351 EZ]- The set of A-regular 
functions with singularities of the above form is closed under addition, multiplication, 
differentiation, and integration |30j . 

Hi) The case of a finite number of singularities on the circle of convergence can be treated 
by a straightforward extension of the above result [331 13"T] . 

Lemma Q] implies a particular singular behaviour of the factorial moment generating 
functions, if the numbers p mj „ satisfy certain typical asymptotic estimates. We write 
(a)k — a ■ (a — 1) • . . . • (a — k + 1) to denote the lower factorial. 

Proposition 2. For m,n G No, let real numbers p m ,n be given. Assume that the numbers 
Vm,n asymptotically satisfy, for k G No, 

— ^2(n) kPmtn ~ A k x- m m lk - 1 (m -> 00), (5) 

n 

for real numbers A k) x c ^ k , where A k > 0, x c > 0, and ^ k {0, —1, —2, . . .}. 
Then, the factorial moment generating functions g k (x) satisfy 

where f k = A k T(j k ). □ 
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Model 
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bargraph polygons 










column-convex polygons 
directed column-convex polygons 


2 
3 


1 

3 


i 

2 


Airy 


diagonally convex directed polygons 










rooted self-avoiding polygons* 










directed convex polygons 


2 
3 


1 

3 


1 

2 


meander 


diagonally convex polygons* 






1 

2 




three-choice polygons 












Table 1: Exponents and area limit laws for prominent polygon models. An asterisk denotes a 
numerical analysis. 



Remarks, i) The above assumption on the growth of the coefficients in Eq. ([5]) is typical 
for polygon models, with 7^ = {k — 6) /(ft, and (f> > 0. 

ii) If the numbers p m>n satisfy, in addition to Eq.([5]), Condition i) of Assumption \T\ this 
implies for exponents of the form 7^ = (k — 9)/<f), where <fi > 0, the estimate 1/2 < < 1. 
raj The proposition implies that the singular part of the factorial moment generating 
function gk(x) is asymptotically equal to the singular part of the corresponding (ordinary) 
moment generating function, 



dq 



(sing) 




P{x,q) 



{sing) 



We give a list of exponents and area limit distributions for a number of polygon models. 
An asterisk denotes that corresponding results rely on a numerical analysis. It appears that 
the value (8, <fi) = (1/3, 2/3) arises for a large number of models. Furthermore, the exponent 
70 seems to determine the area limit law. These two observations will be explained in the 
following section. 



3 Limit distributions 

In this section, we will concentrate on models of square lattice polygons in the fixed perime- 
ter ensemble, and analyse their area law. The uniform ensemble is of particular interest, 
since non-Gaussian limit laws usually appear, due to expected phase transitions at q — 1. 



9 



For non- uniform ensembles q ^ 1, Gaussian limit laws are expected, due to the absence of 
phase transitions. 

There are effective techniques for the uniform ensemble, since the relevant generating 
functions are typically algebraic. This is different from the fixed area ensemble, where 
singularities are more difficult to analyse. It will turn out that the dominant singularity of 
the perimeter generating function determines the limiting area law of the model. We will 
first discuss several examples with different type of singularity. Then, we will describe a 
general result, by analysing classes of g-difference equations (sec e.g. [103]), which exactly 
solvable polygon models obey. Whereas in the case q ^ 1 their theory is developed to some 
extent, the case q = 1 is more difficult to analyse. Motivated by the typical behaviour of 
polygon models, we assume that a g-difference equation reduces to an algebraic equation 
as q approaches unity, and then analyse the behaviour of its solution about q = 1. 

Useful background concerning a probabilistic analysis of counting parameters of combi- 
natorial structures can be found in [23 Ch IX] . See [SU1 Ch 1] and 0, Ch 1] for background 
about asymptotic expansions. For properties of formal power series, see (33 Ch 1.1]. A 
useful reference on the Laplace transform, which will appear below, is |23j . 

3.1 An illustrative example: Rectangles 
3.1.1 Limit law of area 

Let p m ,n denote the number of rectangles of half-perimeter m and area n. Consider the 
uniform fixed perimeter ensemble, with a discrete random variable of area X m defined by 

P(X m = n)= (7) 
The fc-th moments of X m are given explicitly by 

m— 1 



E[X*] = X)(J(m-Q)*— 



m — 1 

i=i 

~m 2k f (x(l - x)) k dx = — t— ^-jTT m 2k (m -> oo), 
Jo (2k + 1)' 

where we approximated the Riemann sum by an integral, using the Euler-MacLaurin sum- 
mation formula. Thus, the random variable X m has mean \i m ~ m 2 /6 and variance 
cr^ ~ m 4 /180. Since the sequence of random variables (X m ) does not satisfy the con- 
centration property linim^oo a m //i m = 0, we expect a non-trivial limiting distribution. 
Consider the normalised random variable 

X m = ^ = 4% (8) 

Since the moments of X m converge as m — » oo, and the limit sequence M k := lmv^oo E[X^] 
satisfies the Carleman condition ^2 k (M 2 k)~ 1 ^ 2k ^ = oo, they define [T7J Ch 4.5] a unique 
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random variable X with moments M^. Its moment generating function M(t) = E[e tx ] is 
readily obtained as 



k=0 v 



The corresponding probability distribution p(x) is obtained by an inverse Laplace trans- 
form, and is given by 

< x < 1 



*>(*) = (V (») 

This distribution is known as the beta distribution /3i,i/2- Together with [T71 Thm 4.5.5], 
we arrive at the following result. 

Theorem 1. The area random variable X m of rectangles Eq. |?p has mean \i m ~ m 2 /Q 
and variance ~ m 4 /180. T7ie normalised random variables X m Eq. converge in 
distribution to a continuous random variable with limit law Pi,i/2 Eq. (Tj|). We also have 
moment convergence. □ 

3.1.2 Limit law via generating functions 

We now extract the limit distribution using generating functions. Whereas the derivation 
is less direct than the previous approach, the method applies to a number of other cases, 
where a direct approach fails. Consider the half-perimeter and area generating function 
P(x, q) for rectangles, 

P(x,q) = ^p m , n x m q n - 

m,n 

The factorial moments of the area random variable X m Eq. (j7]) are obtained from the 
generating function via 

gk 



where (a)k = a ■ (a — 1) • . . . • (a — k + 1) is the lower factorial. The generating function 
P(x,q) satisfies [871 Eq. 5.1] the linear g-difference equation [103] 

P(x, q) = x 2 qP(qx, q) + + ^ . (10) 

1 — qx 

Due to the particular structure of the functional equation, the area moment generating 
functions 

1 d k 

9k ^ x) = k\ W Pi - X ' q) 



9=1 
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are rational functions and can be computed recursively from the functional equation, by 
repeated differentiation w.r.t. q and then setting q—1. (Such calculations are easily per- 
formed with a computer algebra system.) This gives, in particular, 



x 2 , N x 2 



9o(x) = T , rg, 9i(x) 



(1-x) 2 ' yiV ' (1-x) 4 ' 

2x 3 , . 6a; 4 



92{x) = — , g 3 {x 



l-x) b " v ' (1-x) 



; ' 



, , x 4 (l + 22x + x 2 ) , , 12x 5 (l + 82 + x 2 

9A X ) = — r, vin ' 9b{x) 



(1-x) 10 ' ^ ' (1-x) 12 

Whereas the exact expressions get messy for increasing k, their asymptotic form about 
their singularity x c = 1 is simply given by 

9*( x )~ (1 _ x ) 2fc+ 2 (^^1)- (11) 

The above result can be inferred from the functional equation, which induces a recursion 
for the functions gk{x), which in turn can be asymptotically analysed. This method is 
called moment pumping [36]. Below, we will extract the above asymptotic behaviour by 
the method of dominant balance. 

The asymptotic behaviour of the moments of X m can be obtained from singularity 
analysis of generating functions, as described in Lemma [21 Using the functional equation, 
it can be shown that all functions gk{x) are Laurent series about x = 1, with a finite number 
of terms. Hence the remark following Lemma [2] implies for the (factorial) moments of the 
random variable X m Eq. (|S]) the expression 

E[(X m ) k ] E[(X m ) k ] k\ k\ 

~ — — — - ~ = I m, — > oo I 

k\ k\ T{2k + 2) {2k+l)\ 1 h 

in accordance with the previous derivation. 

On the level of the moment generating function, an application of Watson's lemma [HI 
Sec 4.1] shows that the coefficients k\ in Eq. ( flTi) appear in the asymptotic expansion of a 
certain Laplace transform of the (entire) moment generating function E[e~* x ], 



fc>0 

Note that the r.h.s. is formally obtained by term-by-term integration of the l.h.s.. 

Using the arguments of [161 Ch 8.11], one concludes that there exists an s > 0, such 
that there is a unique function F(s) analytic for 3?(s) > s with the above asymptotic 
expansion. It is given by 

F(s) = Ei(s 2 )e s \ (12) 
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where Ei(z) = e tz /tdt is the exponential integral. The moment generating function 
M(£) = E[e~* x ] of the random variable X is given by an inverse Laplace transform of F(s), 

POO 

/ e~ st M(t 2 )tdt = F(s). 
Jo 

Since there are effective methods for computing inverse Laplace transforms [23], the 
question arises whether the function F(s) can be easily obtained. It turns out that the 
functional equation Eq. (|10p induces a differential equation for F(s). This equation can be 
obtained in a mechanical way, using the method of dominant balance. 



3.1.3 Dominant balance 

For a given functional equation, the method of dominant balance consists of a certain 
rescaling of the variables, such that the quantity of interest appears in the expansion of a 
rescaled variable to leading order. The method was originally used as an heuristic tool in 
order to extract the scaling function of a polygon model [M] (see the following section). In 
the present framework, it is a rigorous method. 

Consider the half-perimeter and area generating function P(x, q) as a formal power 
series. The substitution q = 1 — e' is valid, since the coefficients of the power series P(x, q) 
in x are polynomials in q. We get the power series in e, 



k>0 



whose coefficients (— l) k Qk(x) are power series in x. The functional equation Eq. ffTU]) 
induces an equation for H(x,e), from which the factorial area moment generating functions 
Qk{x) may be computed recursively. 

Now replace gk{x) by its expansion about x — 1, 

\ -* fk,i 



9k{x) 



(I — x \2k+2~l ' 
l>0 V ' 



Introducing s — 1 — x, this leads to a power series E(s,Z) in e, 

k>0 \l>0 J 

whose coefficients are Laurent series in H. As above, the functional equation induces an 
equation for the power series E(s,e) in e, from which the expansion coefficients may be 
computed recursively. 

We infer from the previous equation that 

£ (^ 2 ) = (B- 1 )*^) «■ = ^(.,«). (13) 

l>0 \k>0 J 
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Write F(s, e) = J2i>o ^i( s ) el - By construction, the (formal) series F (s) = F(s, 0) coincides 
with the asymptotic expansion of the desired function F(s) Eq. f|T2|) about infinity. 

The above example suggests a technique for computing F (s). The functional equation 
Eq. ( jTUl) for P(x, q) induces, after reparametrisation, differential equations for the functions 
Fi(s), from which Fq(s) may be obtained explicitly. These may be computed by first writing 

P(l ,,) = _i_ F (^LL^, (!_,)■/»), (14) 

and then introducing variables s and e, by setting x = 1 — se and q = 1 — e 2 . Expand the 
equation to leading order in e. This yields, to order e°, the first order differential equation 

si^(s) + 2 - 2s 2 F (s) = 0. 

The above equation translates into a recursion for the coefficients fk,o, from which fkfl = 
can be deduced. In addition, the equation has a unique solution with the prescribed 
asymptotic behaviour Eqn. ( {TBI , which is given by Fq(s) = Ei(s 2 ) e s . 

As we will argue in the next section, Eq. (fT4j) is sometimes referred to as a scaling 
Ansatz, the function F(s,0) appears as a scaling function, the functions Fi(s), for / > 1, 
appear as correction-to-scaling functions. In our formal framework, where the series Fi(s) 
are rescaled generating functions for the coefficients fk,i, their derivation is rigorous. 



3.2 A general method 

In the preceding two subsections, we described a method for obtaining limit laws of counting 
parameters, via a generating function approach. Since this method will be important in 
the remainder of this section, we summarise it here. Its first ingredient is based on the 
so-called method of moments [TTl Thm 4.5.5]. 

Proposition 3. Form,n E No, let real numbers p m ^ n be given. Assume that the numbers 
Pm,n asymptotically satisfy, for k G N 0j 

77 ^{n)kPm,n ~ Ax^m 7 ^ 1 (m -> oo), (15) 

n 

where A^ are positive numbers, and 7fc = (k — 0)/<f>, with real constants 9 and <fi > 0. 
Assume that the numbers Mk := Ak/A satisfy the Carleman condition 

oo 

$>f 2fc )- 1/(2fc) = +oo. (16) 

k=l 

Then the following conclusions hold. 
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i) For almost all m, the random variables X„ 



nX m = n) = (17) 

2 'n Vm,n 



are well defined. We have 



m i/<P 



for a unique random variable X with moments M k , where d denotes convergence in 
distribution. We also have moment convergence. 

ii) If the numbers M/~ satisfy for all t EM. the estimate 

M k t k , , 

lim -±- = 0, (19) 

fe^oo k\ 

then the moment generating function M(t) = K[e~ tx ] of X is an entire function. 
The coefficients AkT(^k) are related to M(t) by a Laplace transform which has, for 
9 > 0, the asymptotic expansion 



oo 

e~ st 



o 



,fc>0 / 



^(-i) fc A fe r( 7fc ) s -^ (s^oc) 



(20) 



u fc>0 



sketch. A straightforward calculation using Eq. ffl5l) leads to 



n(x m ) k ] a_ 

k\ ~ A 



k m kl4 > (m -> oo) 



This implies that the same asymptotic form holds for the (ordinary) moments E[(X m ) fc ]. 
Due to the growth condition Eq. (TIBI) . the sequence (M&) defines a unique random variable 
X with moments M k . Also, moment convergence of the sequence (X m ) to X implies 
convergence in distribution, see [TTJ Thm 4.5.5]. Due to the growth condition Eq. ffl9l) . the 
function M(t) is entire. Hence the conditions of Watson's Lemma [5l Sec 4.1] are satisfied, 
and we obtain Eq. (j20]L □ □ 

Remarks, i) The growth condition Eq. fTl9l) implies the Carleman condition Eq. (fT6]) . All 
examples below have entire moment generating functions M(t). 

ii) If 70 < 0, a modified version of Eq. fl2"01 can be given, see for example staircase polygons 
below. 

Proposition [2] states that assumption Eq. ([15]) translates, at the level of the half- 
perimeter and area generating function P(x,q) = J2 mn Pm, n x m q n , to a certain asymptotic 
expression for the factorial moment generating functions 



1 d k 

9k{x) = y.dq^ P{X,q) 
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Their asymptotic behaviour follows from Eq. (|T5|) . and is 




fk 



(x / x c ) 



(1 - x/x c )^ 



where fk = AkT^k)- Adopting the generating function viewpoint, the amplitudes fk 
determine the numbers Ak, hence the moments Mk = Ak/A Q of the limit distribution. The 
series F(s) = J2 k >o(~ ^-) h fk s ~ lk will be of central importance in the sequel. 

Definition 1 (Area amplitude series). Let Assumption U\ be satisfied. Assume that the 
generating function P(x,q) = ^2 mn Pm,nX m q n satisfies asymptotically 



is called the amplitude series. 

Remarks, i) Proposition [3] states that the area amplitude series appears in the asymptotic 
expansion about infinity of a Laplace transform of the moment generating function of the 
area limit distribution. The probability distribution of the limiting area distribution is 
related to F(s) by a double Laplace transform. 

ii) For typical polygon models, all derivatives of P(x,q) w.r.t. q, evaluated at q = 1, exist 
and have the same radius of convergence, see Proposition [TJ Typical polygon models do 
have factorial moment generating functions of the above form, see the examples below. 

The second ingredient of the method consists in applying the method of dominant 
balance. As described above, this may result in a differential equation (or in a difference 
equation [90J ) for the function F(s). Its applicability has to be tested for each given type 
of functional equation. Typically, it can be applied if the factorial area moment generating 
functions gk(x) Eq. ([!]) have, for values x < x c , a local expansion about x = x c of the form 



where jk,i — (k — 0i)/4> an d > &i- If a transfer theorem such as Lemma [2] applies, 
then the differential equation for F(s) induces a recurrence for the moments of the limit 
distribution. If the differential equation can be solved in closed form, inverse Laplace 
transform techniques may be applied in order to obtain explicit expressions for the moment 
generating function and the probability density. Also, higher order corrections to the 
limiting behaviour may be analysed, by studying the functions Fi(s), for I > 1. See [ST] 
for examples. 




with exponents jk 4- {0> — 1> — 2, . . .}. Then, the formal series 
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3.3 Further examples 



Using the general method as described above, area limit laws for the other exactly solved 
polygon models can be derived. A model with the same area limit law as rectangles 
is convex polygons, compare [87]. We will discuss some classes of polygon models with 
different area limit laws. 

3.3.1 Ferrers diagrams 

In contrast to the previous example, the limit distribution of area of Ferrers diagrams is 
concentrated. 

Proposition 4. The area random variable X m of Ferrers diagrams has mean \i m ~ m 2 /8. 
The normalised random variables X m Eq. iH~8\) converge in distribution to a random variable 
with density p(x) = 8(x — 1/8). 

Remark. It should be noted that the above convergence statement already follows from 
the concentration property linim-^oo a m / ji m = 0, with <r^ ~ m 3 /48 the variance of X m , by 
an explicit analysis of the first three factorial moment generating functions. (By Cheby- 
shev's inequality, the concentration property implies convergence in probability, which in 
turn implies convergence in distribution.) For illustrative purposes, we follow a different 
route via the moment method in the following proof. 



Proof. Ferrers diagrams, counted by half-perimeter and area, satisfy the linear g-difference 
equation [87J Eq (5.4)] 



The perimeter generating function go(x) = x 2 /(l — 2x) is obtained by setting q — 1 in 
the above equation. Hence x c = 1/2. Using the functional equation, it can be shown by 
induction on k that all area moment generating functions gk{x) are rational in go(x) and 
its derivatives. Hence all gk{x) are rational functions. Since the area of a polygon grows 
at most quadratically with the perimeter, we have a bound on the exponent, 7& < 2k + 1, 
of the leading singular part of gk{x). Given this bound, the method of dominant balance 
can be applied. We set 



and introduce new variables s and e by q = 1 — e 2 and 2x = 1 — se. Then an expansion 
of the functional equation yields, to order e°, the ODE of first order F'(s) = 4sF(s) — 1, 
whose unique solution with the prescribed asymptotic behaviour is 
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It can be inferred from the differential equation that all coefficients in the asymptotic 
expansion of F(s) at infinity are nonzero. Hence, the above exponent bound is tight. 
It can be inferred from the functional equation by induction on k that each gk{x) is a 
Laurent polynomial about x c = 1/2. Thus, Lemma [2] applies, and we obtain the moment 
generating function of the corresponding random variable Eq. ( JT8l) as M(s) = exp(— s/8). 
This is readily recognised as the moment generating function of a probability distribution 
concentrated at x — 1/8. □ □ 

A sequence of random variables, which satisfies the concentration property, often leads 
to a Gaussian limit law, after centering and suitable normalisation. This is also the case 
for Ferrers diagrams. 



Theorem 2 ([97]). The area random variable X m of Ferrers diagrams has mean \i m ~ 
m 2 /8 and variance cr^ ~ m 3 /48. The centred and normalised random variables 

X m = X ™~»™ , (21) 

converge in distribution to a Gaussian random variable. □ 

Remarks, i) It is possible to prove this result by the method of dominant balance. The 
idea of proof consists in studying the functional equation of the generating function for the 
"centred coefficients" p m ^ n — /i m . 

ii) The above arguments can also be applied to stack polygons to yield the concentration 
property and a central limit theorem. 

3.3.2 Staircase polygons 

The limit law of area of staircase polygons is the Airy distribution. This distribution (see 
and the survey [52]) is conveniently defined via its moments. 



Definition 2 (Airy distribution [31]). The random variable Y is said to be Airy distributed 
if 

E[Y k ] _ r( 7o ) fc 
k\ r(7 fc )0o' 

where 7& = 3k /2 — 1/2, and the numbers <ftk satisfy, for k > 1, the quadratic recurrence 

1 k 

7fe-i0fe-i + 2 fofa- 1 = °> 

1=0 

with initial condition O = — 1- 

Remarks ( |34L 58J). i) The first moment is E[Y] = \pK. The sequence of moments can 
be shown to satisfy the Carleman condition. Hence the distribution is uniquely determined 
by its moments. 
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ii) The numbers <pk appear in the asymptotic expansion of the logarithmic derivative of 
the Airy function at infinity, 

Al og Ai( S )~^(-l) fc ^ S -^ (a-oo), 

fc>0 

where Ai(x) = ^ J °° cos(t 3 /3 + tx) dt is the Airy function. 

Hi) Explicit expressions for the numbers <pk are known [58]. They are, for k > 1, given by 



oo x 3(k-l)/2 

dx, 



^ K 4vr 2 J Ai(x) 2 + Bi(x) 2 

where Bi(z) is the second standard solution of the Airy differential equation f"(z) — zf(z) = 
0. 

iv) The Airy distribution appears in a variety of contexts [34J. In particular, the random 
variable Y j v8 describes the law of the area of a Brownian excursion. See also [76] for an 
overview from a physical perspective. 

Explicit expressions have been derived for the moment generating function of the Airy 
distribution and for its density. 

Fact 1 ([TjH [661 EH El])- The moment generating function M{t) = E[e~* y ] of the Airy 
distribution satisfies the modified Laplace transform 

j=[ (.- -DM^/V^d* = 2V3 _ . (22) 

T/ie moment generating function M(t) is given explicitly by 

oo 

M(2~ 3 / 2 t) = V2^t ]Texp (-^t 2 /^" 1 / 3 ) , 
fc=i 

where the numbers —[3k are the zeros of the Airy function. Its density p[x) is given explicitly 
by 

k=i ^ ' 

where Vk = 2/5|/(27a; 2 ) and U(a,b,z) is the confluent hypergeometric function. □ 

Remarks, i) The confluent hypergeometric function U(a,b; z) is defined as [1] 

, 7T / i : Fi[a 1 M ^^[l + a-fe^-fe^] 

la ' ,Z) smirb\T(l + a-b)T(b) T(a)V(2-b) 

where i-Fi[a; 6; z] is the hypergeometric function 

„ r , , a z a(a + 1) z 2 

1 i, l[a;M ] = i + -- + ^ TT i- + ... 
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ii) The moment generating function and its density are obtained by two consecutive inverse 
Laplace transforms of Eq. (122]) , see [671 EE] and [991 EI] • 

Hi) In the proof of the following theorem, we will derive Eq. (J22l using the model of staircase 
polygons. This shows, in particular, that the coefficients <pk appear in the asymptotic 
expansion of the Airy function. 

Theorem 3. The normalised area random variables X m of staircase polygons Eq. [W) 
satisfy 

A m d Y . . 
0F/4 ^ 

where Y is Airy distributed according to Definition^ We also have moment convergence. 

Remark. Given the functional equation of the half-perimeter and area generating function 
of staircase polygons, 

2 

P{x > q) = l-2 xq X -P( qX , q ) (23) 

(see [88] for a recent derivation), this result is a special case of Theorem H] below, which is 
stated in [251. 



Proof. We use the method of dominant balance. From the functional equation Eq. (1231) . 
we infer go(x) = 1/4 + a/1 — 4x/2 + (1 — 4x)/4. Hence x c = 1/4. The structure of the 
functional equation implies that all functions gk{x) can be written as Laurent series in 
s = y/1 — Ax, see also Proposition [7J below. Explicitly, we get gi(x) = x 2 /(l — Ax). This 
suggests 7^ = (3k — l)/2. An upper bound of this form on the exponent 7 fc can be derived 
without too much effort from the functional equation, by an application of Faa di Bruno's 
formula, see also [SSI Prop (4.4)]. Thus, the method of dominant balance can be applied. 
We set 

P(x,,) = i + (l-,)V3 F (_i^i|_, (1 - 5) ./3j 

and introduce variables s, e by 4i = 1 — se 2 and q = 1 — e 3 . In the above equation, 
we excluded the constant 1/4 =: P^ re9 '(x,q), since it does not contribute to the moment 
asymptotics. Expanding the functional equation to order e 2 gives the Riccati equation 

F'(s) +AF(s) 2 - s = 0. (24) 

It follows that the coefficients fk of F(s) satisfy, for k > 1, the quadratic recursion 

k 

Tfc-iA-i + 4 ^2 hf*- 1 = °' 

1=0 

with initial condition /o = —1/2. A comparison with the definition of the Airy distribution 
shows that <pk = 2 2k+1 fk- Using the closure properties of A- regular functions, it can be 
inferred from the functional equation that (the analytic continuation of) each factorial 
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moment generating function gt{x) is A-regular, with x c = 1/4, see also Proposition [7J 

below. Hence the transfer theorem Lemma [2] can be applied. We obtain 4X m Y in 
distribution and for moments, where Y is Airy distributed. □ □ 

Remarks, i) The unique solution F(s) of the differential equation in the above proof 
Eq. (1241) . satisfying the prescribed asymptotic behaviour, is given by 

F{ 8 ) = ~]ogM{A 1 ^ 8 ). (25) 

The moment generating function M(t) of the limiting random variable X = lim ni _ >00 X m 
is related to the function F(s) via the modified Laplace transform 

(e~ st -l)M(t 3 / 2 )^dt = A^{F{s) - F(0)), 

where the modification has been introduced in order to ensure a finite integral about the 
origin. This result relates the above proof to Proposition [IJ 

ii) The method of dominant balance can be used to obtain corrections Fi(s) to the limiting 
behaviour [87] . 

The fact that the area law of staircase polygons is, up to normalisation, the same as 
that of the area under a Brownian excursion, suggests that there might be a combinato- 
rial explanation. Indeed, as is well known, there is a bijection [2TJ [98] between staircase 
polygons and Dyck paths, a discrete version of Brownian excursions [2] , see figure [2] [88J . 
Within this bijection, the polygon area corresponds to the sum of peak heights of the Dyck 
path, but not to the area below the Dyck path. For more about this connection, see the 




Figure 2: |88j A combinatorial bijection between staircase polygons and Dyck paths |21l [98] . 
Column heights of a polygon correspond to peak heights of a path. 

remark at the end of the following subsection. 
3.4 g-difference equations 

All polygon models discussed above have an algebraic perimeter generating function. More- 
over, their half-perimeter and area generating function satisfies a functional equation of 
the form 

P(x, q) = G(x, q, P(x, q), P(qx, q)), 
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for a real polynomial G(x,q,yo,yi). Since, under mild assumptions on G, the equation 
reduces to an algebraic equation for P(x, 1) in the limit q — > 1, it may be viewed as a 
"deformation" of an algebraic equation. In this subsection, we will analyse equations of 
this type at the special point (x,q) = (x c , 1), where x c is the radius of convergence of 
P(x, 1). It will appear that the methods used in the above examples also can be applied 
to this more general case. 

The above equation falls into the class of q-difference equations [103J. While particular 
examples appear in combinatorics in a number of places, see e.g. [37J, the asymptotic 
behaviour of equations of the above form seems to have been systematically studied initially 
in [251 187] . The study can be done in some generality, e.g., also for non-polynomial power 
series G, for replacements more general than x i— > qx, and for multivariate generalisations, 
see [89] and [25]. For simplicity, we will concentrate on polynomial G, and then briefly 
discuss generalisations. Our exposition closely follows [891 EZ] • 



3.4.1 Algebraic g-difference equations 

Definition 3 (Algebraic g-difference equation J25j[87]). An algebraic q-difference equation 
is an equation of the form 

P{x, q) = G{x, q, P{x, q), P{qx, g), . . . , P(q N x, g)), (26) 

where G(x, q, yo, y\, . . . , yw) is a complex polynomial. We require that 

dG 

G(0,g,0,0,...,0) = 0, — (0,g, 0,0,..., 0) = (k = 0, 1, . . . , N). 

cyk 

Remarks, i) See [103] for an overview of the theory of g-difference equations. As q 
approaches unity, the above equation reduces to an algebraic equation. 
ii) Asymptotics for solutions of algebraic g-difference equations have been considered in 
[25] . The above definition is a special case of [891 Def 2.4], where a multivariate extension 
is considered, and where G may be non-polynomial. Also, replacements more general than 
x i— > f(q)x are allowed. Such equations are called g-functional equations in [HH]. The 
results presented below apply mutatis mutandis also to g-functional equations. 

The algebraic g-difference equation in Definition [3] uniquely defines a (formal) power 
series P(x, g) satisfying P(0, g) = 0. This is shown by analysing the implied recurrence 
for the coefficients p m {q) of P(x,q) = 'Yj m >oVrn{.q)x m , see also [89, Prop 2.5]. In fact, 
Pm(q) is a polynomial in q. The growth of its degree in m is not larger than cm 2 for some 
positive constant c, hence the counting parameters are rank 2 parameters [25]. In our 
situation, such a bound holds, since the area of a polygon grows at most quadratically 
with its perimeter. 

From the preceding discussion, it follows that the factorial moment generating functions 

1 d k 

9k ( X > = k\dq* P{ - X " q) 
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are well-defined as formal power series. In fact, they can be recursively determined from the 
g-difference equation by implicit differentiation, as a consequence of the following proposi- 
tion. 



Proposition 5 f [87[ [89]). Consider the derivative of order k > of an algebraic q- 
difference equation Eq. ([251) w.r.t. q, evaluated at q = 1. It is linear in gk(x), and its 
r.h.s. is a complex polynomial in the power series gi(x) and its derivatives up to order 



Remarks, i) This statement can be shown by analysing the /c-th derivative of the q- 
difference equation, using Faa di Bruno's formula [18] . 

ii) It follows that every function gu(x) is rational in gi(x) and its derivatives up to order 
k — l, where < I < k. Since G is a polynomial, gk(x) is algebraic, by the closure properties 
of algebraic functions. 

We discuss analytic properties of the (analytic continuations of the) factorial moment 
generating functions gk{x). These are determined by the analytic properties of go(x) = 
P(x,l). We discuss the case of a square-root singularity of P(x, 1), which often occurs 
for combinatorial structures, and which is well studied, see e.g. [79J Thm 10.6] or [37J 
Ch VII. 4]. Other cases may be treated similarly. We make the following assumption: 

Assumption 2. The q-difference equation in Definition^ has the following properties: 

i) All coefficients of the polynomial G(x, q, y , yi, . . . , y^) are non-negative. 

ii) The polynomial Q(x, y) := G(x, 1, y, y, . . . , y) satisfies Q(x, 0) ^ and has degree at 
least two in y. 

Hi) P(x, 1) = J2 m >iPm xm is aperiodic, i.e., there exist indices 1 < i < j < k such that 
PiPjPk 7^ 0, while gcd(j — i, k — i) = 1. 

Remarks, i) The positivity assumption is natural for combinatorial constructions. There 
are, however, g-difference equations with negative coefficients, which arise from systems of 
g-difference equations with non-negative coefficients by reduction. Examples are convex 
polygons [87J Sec 5.4] and directed convex polygons, see below. 

ii) Assumptions i) and ii) result in a square-root singularity as the dominant singularity 



Hi) Assumption Hi) implies that there is only one singularity of P(x, 1) on its circle of 
convergence. Since P(x, 1) has non-negative coefficients only, it occurs on the positive real 
half-line. The periodic case can be treated by a straightforward extension |37j . 

An application of the (complex) implicit function theorem ensures that P(x, 1) is an- 
alytic at the origin. It can be analytically continued, as long as the defining algebraic 
equation remains invertible. Together with the positivity assumption, one can conclude 
that there is a number < x c < oo, such that the analytic continuation of P(x, 1) satisfies 
y c = \im x y Xc P(x, 1) < oo, with 



k — l, where I = 0, . . . , k. 



□ 



of P{x,l). 



Q(x c ,y c ) = y, 



d_ 

dy 



Q(x c ,y) 



l. 
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With the positivity assumption on the coefficients, it follows that 

> 0. (27) 



1 2 



B: =2W Q(Xc,y) 



d 

>0, C := —Q(x,y c ) 
y=yc ( '' 



These conditions characterise the singularity of P(x, 1) at root. It 

can be shown that there exists a locally convergent expansion of P(x, 1) about x = x c , and 
that P(x, 1) is analytic for \x\ < x c . We have the following result. Recall that a function 
f(z) is A-regular if it is analytic in the indented disc A = {z : \z\ < x c + r], \ Aig(z — x c )\ > 
4>} for some r\ > and some 0, where < < ir/2. 

Proposition 6 ([751 EH EH])- Given Assumption IE, the power series P(x,l) is analytic 
at x = 0, with radius of convergence x c . Its analytic continuation is A-regular, with a 
square-root singularity at x = x c and a local Puiseux expansion 



P(x, 1) = y c + ^2 ^o.^ 1 ~ x l x ^ 



1/2+1/2 

j 



1=0 



where y c = \im x y Xc P(x,l) < oo and / ,o — —\fx c C/B, for constants B > and C > 
as in Eq. pT7] ). The numbers foj can be recursively determined from the q- difference 
equation. □ 

The asymptotic behaviour of P(x, 1) = go{x) carries over to the factorial moment 
generating functions gk{x). 

Proposition 7 ([89]). Given Assumption 0, all factorial moment generating functions 
gk{x) are, for k > 1, analytic at x = 0, with radius of convergence x c . Their analytic 
continuations are A-regular, with local Puiseux expansions 



1=0 

where jk = 3k/2 — 1/2. The numbers /^q = fk are, for k > 2, characterised by the 
recursion 

1 k 

lk-lfk-l + -TJ- ^ flfk-l = 0, 
1=0 

and the numbers fo<0 and f\ > are given by 



[Cx~ c Ef=ifcff(^,l,l/ c ,y c ,---,yc) 
/o = "V~B"' h = B ' ( } 

for constants B > and C > as in Eq. pT7|). □ 
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Remarks, i) This result can be obtained by a direct analysis of the g-difference equation, 
applying Faa di Bruno's formula, see also [HTJ Sec 2.2]. 

ii) Alternatively, it can be obtained by applying the method of dominant balance to the 
g-difference equation. To this end, one notes that all functions gk{x) are Laurent series in 
y/l — x/x c , and that their leading exponents are bounded from above by 7^. (An upper 
bound on an exponent is usually easier to obtain than its exact value, since cancellations 
can be ignored). With these two ingredients, the method of dominant balance, as described 
above, can be applied. The differential equation of the function F(s) then translates, via 
a transfer theorem, into the above recursion for the coefficients. See [SHI Sec 5]. 

The above result can be used to infer the limit distribution of area, along the lines of 
Section El 

Theorem 4 ([25J EH])- Let Assumption^ be satisfied. For the solution of an algebraic 
q-difference equation P(x,q) = 'Yl imn Pm,,nX m q n , let X m denote the random variable 

= n) 



Pm,n 

(which is well-defined for almost all m). The mean of X m is given by 

E[X m ] ~ 2v^A m 3/ 2 (774 ► 00) , 
l/o I 

where the numbers /o and f\ are given in Eq. (fD5p. The sequence of normalised random 
variables X m converges in distribution, 

X m d Y 
X m = — = > —= [m — > 00), 

'7T 



E[X n 

where Y is Airy distributed according to Definition^ We also have moment convergence. 

□ 

Remarks, i) An explicit calculation shows that (p^ = \fo\~ 1 (27^) fk- Together with 
Proposition ^\ the claim of the proof follows by standard reasoning, as in the examples 
above. 

ii) The above theorem appears in [251 Thm 3.1], together with an indication of the ar- 
guments of a proof. [There is a misprint in the definition of 7 in [251 Thm 3.1]. In our 
notation 7 = 4Bf\.] Within the more general setup of g-functional equations, the theorem 
is a special case of [891 Thm 1.5]. 

Hi) The above theorem is a kind of central limit theorem for combinatorial constructions, 
since the Airy distribution arises under natural assumptions for a large class of combina- 
torial constructions. For a connection to certain Brownian motion functionals, see below. 
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3.4.2 g-functional equations and other extensions 



We discuss extensions of the above result. Generically, the dominant singularity of P(x, 1) 
is a square-root. The case of a simple pole as dominant singularity, which generalises 
the example of Ferrers diagrams, has been discussed in [87J. Under weak assumptions, 
the resulting limit distribution of area is concentrated. Other singularities can also be 
analysed, as shown in the examples of rectangles above and of directed convex polygons in 
the following subsection. Compare also |90j . 

The case of non-polynomial G can be discussed along the same lines, with certain 
assumptions on the analyticity properties of the series G. In the undeformed case q = 1, 
it is a classical result [371 Ch VII. 3] that the generating function has a square-root as 
dominant singularity, as in the polynomial case. One can then argue along the above 
lines that an Airy distribution emerges as the limit law of the deformation variable [89l 
Thm 1.5]. Such an extension is relevant, since prominent combinatorial models, such as the 
Cayley tree generating function, fall into that class. See also the discussion of self-avoiding 
polygons below. 

The above statements also remain valid for more general classes of replacements x i— > qx, 
e.g., for replacements x i— > f(q)x, where f(q) is analytic for < q < 1, with non-negative 
series coefficients about q = 0. More interestingly, the idea of introducing a g-deformation 
may be iterated [25], leading to equations such as 

P(x, qx, . . . , q M ) = G(x, P(xqi ■ . . . ■ q M , <?i<?2 • • • • • qM, <?2<?3 • • • • • Qm, ■ ■ ■ , (29) 

The counting parameters corresponding to q^ are rank k + 1 parameters, and limit distri- 
butions for such quantities have been derived for some types of singularities [77J [78j [88] . 
There is a central limit result for the generic case of a square- root singularity [89]. This 
generalisation applies to counting parameters, which decompose linearly under a combina- 
torial construction. These results can also be obtained by an alternative method, which 
generalises to non-linear parameters, see [oT| . 

The case where the limit q to unity in a g-difference equation is not algebraic, has not 
been discussed. For example, if G(x,q, P(x,q), P(qx,q)) = for some polynomial G, the 
limit q to unity might lead to an algebraic differential equation for P(x, 1). This may be 
seen by noting that 



for f(x) differentiable at x. Such equations are possibly related to polygon models such as 
three-choice polygons [H] or punctured staircase polygons |45j . Their perimeter generating 
function is not algebraic, hence the models do not satisfy an algebraic g-difference equation 
as in Definition [3l 

3.4.3 A stochastic connection 

Lastly, we indicate a link to Brownian motion, which appears in [99, 100J and was further 
developed in [TTJ HH [SH [88] . As we saw in Section 13.21 limit distributions can, under 
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certain conditions, be characterised by a certain Laplace transform of their moment gen- 
erating functions. This approach, which arises naturally from the viewpoint of generating 
functions, can be applied to discrete versions of Brownian motion, excursions, bridges or 
meanders. Asymptotic results are results for the corresponding stochastic objects. In fact, 
distributions of some functionals of Brownian motion have apparently first been obtained 
using this approach [99], HOOj . 

Interestingly, a similar characterisation appears in stochastics for functionals of Brow- 
nian motion, via the Feynman-Kac formula. For example, Louchard's formula [66] relates 
the logarithmic derivate of the Airy function to a certain Laplace transform of the moment 
generating function of the law of the Brownian excursion area. Distributions of functionals 
of Brownian motion can also be obtained by a path integral approach, see [75] for a recent 
overview. 

The discrete approach provides an alternative method for obtaining information about 
distributions of certain functionals of Brownian motion. For such functionals, it provides 
an alternative proof of Louchard's formula [7_7J EE]- It leads, via the method of dominant 
balance, quite directly to moment recurrences for the underlying distribution. These have 
been studied in the case of rank k parameters for discrete models of Brownian motion. 
In particular, they characterise the distributions of integrals over (k — l)-th powers of the 
corresponding stochastic objects (TTj [78j [EH EE]- Such results have apparently not been 
previously derived using stochastic methods. The generating function approach can also be 
applied to classes of g-functional equations with singularities different from those connected 
to Brownian motion. For a related generalisation, see [TO] . 

Vice versa, results and techniques from stochastics can be (and have been) analysed in 
order to study asymptotic properties of polygons. An example is the contour process of 
simply generated trees [38] . which asymptotically describes the area of a staircase polygon. 

see also jsacna HUES]. 

3.5 Directed convex polygons 

We show that the limit law of area of directed convex polygons in the uniform fixed 
perimeter ensemble is that of the area of the Brownian meander. 

Fact 2 ( |100[ Thm 2]). The random variable Z of area of the Brownian meander is char- 
acterised by 

E[Z k ] _ r(q ) uj k 1 
k\ r(a fc ) w 2^' 
•where a k = 3k/2 + 1/2. The numbers u k satisfy for k > 1 the quadratic recurrence 

k 

a k -iuJk-i + ^2 ^~ l<Jt} k-i = 0, 

1=0 

with initial condition luq — 1, where the numbers <f) k appear in the Airy distribution as in 
Definition d □ 
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Remarks, i) This result has been derived using a discrete meander, whose length and 
area generating function is described by a system of two algebraic g-difference equations, 
see (73, Prop 1]. 

ii) We have K[Z] = 3\^2n/8 for the mean of Z. The random variable Z is uniquely 



determined by its moments. The numbers Uk appear in the asymptotic expansion [1001 
Thm 3] 

1 — 3 rf Ai(t) dt v-^ , 

^ ' fc>0 

where Ai(x) = ^ Jq 00 cos(t 3 /3 + tx) dt is the Airy function. 

Explicit expressions have been derived for the moment generating function and for the 
distribution function of Z. 

Fact 3 ([100, Thm 5]). The moment generating function M(t) = K[e~ tz ] of Z satisfies the 
Laplace transform 

poo 1 

/ e- st M(V2t 3/2 )-^dt = 0Fft(s). (30) 

It is explicitly given by 

oo 

M(t) = 2- 1 / 6 t 1 / 3 ^ J R fc exp(-/? fc t 2 / 3 2- 1 /3) 

k=l 

for dt(t) > 0, where the numbers —f3k are the zeroes of the Airy function, and where 

/5 fc (l + 3/ /3fc Ai(-t)dt) 



3Ai'(-&) 



The random variable Z has a continuous density p(y), with distribution function R(x) 
JoP(v) d V 9 wen b V 



oo 

7T 



R{x) = 7W7^E^ e ""% 1/3Ai ((3^/2) 2/3 ^ 



k=l 

where v k = (/3 fc ) 3 /(27x 2 ). □ 

Remark. The moment generating function and the distribution function are obtained by 
two consecutive inverse Laplace transforms of Eq. (1301 . 

Theorem 5. The normalised area random variables X m of directed convex polygons Eq. ( TJ7| ) 
satisfy 

X m — > -Z (m -> oo), 

where Z is the area random variable of the Brownian meander as in Fact\^ We also have 
moment convergence. 
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Proof. A system of g-difference equations for the generating function Q(x, y, q) of directed 
convex polygons, counted by width, height and area, has been given in [9j Lemma 1.1]. It 
can be reduced to a single equation, 

q(qx - l)Q(x, y, q) + ((1 + q)(P(x, y, q)+y)) Q(qx } y, q) + 
(xyq - y 2 + P(x, y, q) (qx - y - 1)) Q(q 2 x, y, q) (31) 
- q 2 xy (y + P(x, y, q) - 1) = 0, 

where P(x, y, q) is the width, height and area generating function of staircase polygons. 
Setting q=l and x = y yields the half-perimeter generating function 

Hence x c = 1/4 for the radius of convergence of Q(x,x, 1). 

It is possible to derive from Eq. (13T1) a g-difference equation for the (isotropic) half- 
perimeter and area generating function Q(x,q) = Q(x,x,q) of directed convex polygons. 
This is due to the symmetry Q(x,y,q) = Q(y,x,q), which results from invariance of the 
set of directed convex polygons under reflection along the negative diagonal y = —x. Since 
this equation is quite long, we do not give it here. By arguments analogous to those of the 
previous subsection, it can be deduced from this equation that all area moment generating 
functions gk(x) of Q(x, 1) are Laurent series in s = y/l — 4x, see also [89, Prop (4.3)]. The 
leading singular exponent of gk(x), defined by gk(x) ~ h k (l — x/x c )~ ak as x f x c , can be 
bounded from above by a k < 3k/2 + 1/2, see also [HHl Prop (4.4)] for the argument. We 
apply the method of dominant balance, in order to prove that a k = 3&/2 + 1/2 and to 
yield recurrences for the coefficients hk- We define 

P(»,,) = J + (i- 8 )v» F (^J!^ Sl (i_ a )V.) i 

Q(x, ,) = (!- <,)-"* H ( J-^lj, (1 - «) 1/3 ) , 

where F(s) = F(s,0) has already been determined in Eq. ( |25l) . We set 4x = 1 — se 2 , 
q = 1 — e 3 , and expand the g-difference equation to leading order in e. We get for H(s) := 
H(s, 0) the inhomogeneous linear differential equation of first order 

H\s)+4H(s)F(s) + l = 0. 

o 

This implies for the coefficients hk of H(s) = ^2 k >o hkS~ ak and fk of F(s) = J2k>o fkS~ lk 
for k > 1 the quadratic recursion 

k 

ak-ihk-! + 4 ^ fi h k-i = °) 

1=0 
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where ho = 1/16. Using fj. = 2~ 2k ~ 1 <pk, we obtain the meander recursion in Fact [2] by 
setting hk = 2~ k ~ 4 ujj t . It can be inferred from the functional equation that (the analytic 
continuations of) all factorial moment generating functions are A-regular, with x c = 1/4. 

Thus Lemma [2] applies, and we conclude X m —> Z/2. □ □ 

Remarks, i) The above theorem states that the limit distribution of area of directed 
convex polygons coincides, up to normalisation, with the area distribution of the Brownian 
meander [100] . This suggests that there might exist a combinatorial bijection to discrete 
meanders, in analogy to that between staircase polygons and Dyck paths. Up to now, 
a "nice" bijection has not been found, see however [51 [72j for combinatorial bijections to 
discrete bridges. 

ii) The above proof relies on a g-difference equation for the isotropic generating function 
Q(x,x,q). Up to normalisation, the meander distribution also appears for the anisotropic 
model Q(x,y,q), where 0<y<l/2is fixed, as can be shown by a considerably simpler 
calculation. The normalisation constant coincides with that of the isotropic model for 
y = 1/2. The latter statement is also a consequence of the fact that the height random 
variable of directed polygons is asymptotically Gaussian, after centering and normalisation. 
Analogous considerations apply to the relation between isotropic and anisotropic versions 
of the other polygon classes. 



3.6 Limit laws away from (x c , 1) 

As motivated in the introduction, limit laws in the fixed perimeter ensemble for q ^ 1 are 
expected to be Gaussian. The same remark holds for the fixed area ensemble for x ^ x c . 
There are partial results for the model of staircase polygons. The fixed area ensemble can, 
for x < x c and q near unity, be analysed using Fact [7] of the following section. For staircase 
polygons in the uniform fixed area ensemble x — 1, the following result holds. 

Fact 4 ([371 Prop IX.ll]). Consider the perimeter random variable of staircase polygons 
in the uniform fixed area ensemble, 

Pfc = m) = 

l^im Pm,n 

The variable Y n has mean \x n ~ fi-n and standard deviation o~ n ~ cry/n, where the numbers 
fM and a satisfy 

[i = 0.8417620156 . . . , a = 0.4242065326 . . . 

The centred and normalised random variables 

in i 
On 

converge in distribution to a Gaussian random variable. □ 
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Remark. The above result is proved using an explicit expression for the half-perimeter 
and area generating function, as a ratio of two g-Bessel functions. It can be shown that this 
expression is meromorphic about (x,q) = (1, q c ) with a simple pole, where q c is the radius 
of convergence of the generating function P(l,q). The explicit form of the singularity 
about (l,<7c) yields a Gaussian limit law. 

There are a number of results for classes of column-convex polygons in the uniform 
fixed area ensemble, typically leading to Gaussian limit laws. The upper and lower shape 
of a polygon can be described by Brownian motions. See [69, [TOj, [Tl] for details. It would 
be interesting to prove convergence to a Gaussian limit law within a more general frame- 
work, such as g-difference equations. Analogous questions for other functional equations, 
describing counting parameters such as horizontal width, have been studied in [23] . 

3.7 Self-avoiding polygons 

A numerical analysis of self-avoiding polygons, using data from exact enumeration [9T1 
92J, supports the conjecture that the limit law of area is, up to normalisation, the Airy 
distribution. 

Let p m ,n denote the number of square lattice self-avoiding polygons of half-perimeter 
m and area n. Exact enumeration techniques have been applied to obtain the numbers 
Pm,n for all values of n for given m < 50. Numerical extrapolation techniques yield very 
accurate estimates of the asymptotic behaviour of the coefficients of the factorial moment 
generating functions. To leading order, these are given by 

[x m ]g k (x) = i ]T(n) fc p m , n ~ A^-m 3 ^ 3 / 2 - 1 (m - oo), (32) 

n 

for positive amplitudes A^. The above form has been numerically checked [91] [92] for values 
k < 10 and is conjectured to hold for arbitrary k. The value x c is the radius of convergence 
of the half-perimeter generating function of self-avoiding polygons. The amplitudes 
have been extrapolated to at least five significant digits. In particular, we have 

x c = 0.14368062927(2), A = 0.09940174(4), A x = 0.0397886(1), 

where the numbers in brackets denote the uncertainty in the last digit. An exact value of 
the amplitude A\ = l/(87r) has been predicted [U] using field-theoretic arguments. 

The particular form of the exponent implies that the model of rooted self-avoiding poly- 
gons p m ,n = mpm,n has the same exponents = 2/3 and 6 = 1/3 as staircase polygons. In 
particular, it implies a square-root as dominant singularity of the half-perimeter generat- 
ing function. Together with the above result for g-functional equations, this suggests that 
(rooted) self-avoiding polygons might obey the Airy distribution as a limit law of area. 

A natural method to test this conjecture consists in analysing ratios of moments, such 
that a normalisation constant is eliminated. Such ratios are also called universal amplitude 
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ratios. If the conjecture were true, we would have asymptotically 

E[x m ]* r( 7t )r( 7o )*-i # 

for the area random variables X m as in Eq. (fl71) . The numbers 0^ and exponents jk are 
those of the Airy distribution as in Definition [2j The above form was numerically confirmed 
for values of k < 10 to a high level of numerical accuracy. The normalisation constant is 
obtained by noting that E[Y] = ^pn. 

Conjecture 1 (cf [9T| 152]). Let p m)n denote the number of square lattice self-avoiding 
polygons of half-perimeter m and area n. Let X m denote the random variable of area in 
the uniform fixed perimeter ensemble, 



F(X m = n 

We conjecture that 



Pr, 



Y2nP^ 
X m d Y 



E[X m ] 

where Y is Airy distributed according to Definition^ 
Remarks, i) Field theoretic arguments [15] yield A\ = 1/(8-71"). 

ii) References [9lT 192] contain conjectures for the scaling function of self-avoiding polygons 
and rooted self-avoiding polygons, see the following section. In fact, the numerical analysis 
in [9Tj [92] mainly concerns the area amplitudes A^, which determine the limit distribution 
of area. 

Hi) The area law of self- avoiding polygons has also been studied [9TJ [92] on the triangular 
and hexagonal lattices. As for the square lattice, the area limit law appears to be the Airy 
distribution, up to normalisation. 

iv) It is an open question whether there are non-trivial counting parameters other than the 
area, whose limit law (in the fixed perimeter ensembles) coincides between self-avoiding 
polygons and staircase polygons. See |88j for a negative example. This indicates that 
underlying stochastic processes must be quite different. 

v) A proof of the above conjecture is an outstanding open problem. It would be interest- 
ing to analyse the emergence of the Airy distribution using stochastic Loewner evolution 
[60] . Self-avoiding polygons at criticality are conjectured to describe the hull of critical 
percolation clusters and the outer boundary of two-dimensional Brownian motion [60J. 

A numerical analysis of the fixed area ensemble along the above lines again shows 
behaviour similar to that of staircase polygons. This supports the following conjecture. 

Conjecture 2. Consider the perimeter random variable of self-avoiding polygons in the 
uniform fixed area ensemble, 

P(n = m)= 

Pra.n 
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The random variable Y n is conjectured to have mean fi n ~ // • n and standard deviation 
cr„ ~ Oyfn, where the numbers \i and o satisfy 

li = 1.855217(1), a 2 = 0.3259(1), 

where the number in brackets denotes the uncertainty in the last digit. The centred and 
normalised random variables _ 

are conjectured to converge in distribution to a Gaussian random variable. 

The above conjectures, together with the results of the previous subsection, also raise 
the question whether rooted square-lattice self-avoiding polygons, counted by half-perimeter 
and area, might satisfy a g-functional equation. In particular, it would be interesting to 
consider whether rooted self-avoiding polygons might satisfy 

P{x) = G(x,P(x)), (33) 

for some power series G(x,y) in x, y. If the perimeter generating function P(x) is not 
algebraic, this excludes polynomials G(x, y) in x and y. Note that the anisotropic perimeter 
generating function of self-avoiding polygons is not D-fmite [EE]- It is thus unlikely that 
the isotropic perimeter generating function is .D-finite and, in particular, algebraic. On the 
other hand, solutions of Eq. (|33|) need not be algebraic nor D-finite. An example is the 
Cayley tree generating function T(x) satisfying T(x) = xexp(T(x)), see [33J. 



3.8 Punctured polygons 

Punctured polygons are self-avoiding polygons with internal holes, which are also self- 
avoiding polygons. The polygons are also mutually avoiding. The perimeter of a punctured 
polygon is the sum of the lengths of its boundary curves, the area of a punctured polygon is 
the area of the outer polygon minus the area of the holes. Apart from intrinsic combinato- 
rial interest, models of punctured polygons may be viewed as arising from two-dimensional 
sections of three-dimensional self-avoiding vesicles. Counted by area, they may serve as an 
approximation to the polyomino model. 

We consider, for a given subclass of self-avoiding polygons, punctured polygons with 
holes from the same subclass. The case of a bounded number of punctures of bounded 
size can be analysed in some generality. The case of a bounded number of punctures of 
unbounded size leads to simple results if the critical perimeter generating function of the 
model without punctures is finite. 

For a given subclass of self-avoiding polygons, the number p m>n denotes the number 
of polygons with half-perimeter m and area n. Let pm,n denote the number of polygons 
with r > 1 punctures whose half-perimeter sum equals s. Let pm,n denote the number of 
polygons with r > 1 punctures of arbitrary size. 
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Theorem 6 ( [Ml Thms 1,2]). Assume that, for a class of self- avoiding polygons with- 
out punctures, the area moment coefficients p^m = ^2 n >o nk Pm,n have, for k G N 0; the 
asymptotic form 

pJ} J ~ A k x- m m? k - 1 (m -> oo), 

for numbers A k > 0, for < x c < 1 and /or 7^ = (fc — #)/</>, where < <p < 1. Let 
go{x) = 'Yli m >oPmX m denote the half-perimeter generating function. 

Then, the area moment coefficient pm^'^ = n k pm$ of the polygon class with r > 1 
punctures whose half-perimeter sum equals s is, for k E No, asymptotically given by 

p^ Ks) ~ A^x^m^- 1 (m -»• 00), 

w/jere = ^f L a;^[x s ](oo(x)) r . 

If6>0, the area moment coefficient Pm^ = Yl n nh P^ n °f the polygon class with r > 1 
punctures of arbitrary size satisfies, for k e No, asymptotically 

p^ k) ~ 4 r)a; c " mm7fc+r_1 (m -> 00), 
where the amplitudes A k are given by 




A k +r(,go{x c )) r 



Remarks. ^ The basic argument in the proof of the preceding result involves an estimate 
of interactions of hole polygons with one another or with the boundary of the external 
polygon, which are shown to be asymptotically irrelevant. This argument also applies in 
higher dimensions, as long as the exponent satisfies < <p < 1. 

ii) In the case of an infinite critical perimeter generating function, such as for subclasses of 
convex polygons, boundary effects are asymptotically relevant, if punctures of unbounded 
size are considered. The case of an unbounded number of punctures, which approximates 
the polyomino problem, is unsolved. 

Hi) The above result leads to new area limit distributions. For rectangles with r punctures 
of bounded size, we get Pr+1,1/2 as the limit distribution of area. For staircase polygons with 
punctures, we obtain generalisations of the Airy distribution, which are discussed in |94j . 
In contrast, for Ferrers diagrams with punctures of bounded size, the limit distribution of 
area stays concentrated. 

iv) The theorem also applies to models of punctured polygons, which do not satisfy an 
algebraic g-difference equation. An example is given by staircase polygons with a staircase 
hole of unbounded size, whose perimeter generating function is not algebraic [43] . 

3.9 Models in three dimensions 

There are very few results for models in higher dimensions, notably for models on the cubic 
lattice. There are a number of natural counting parameters for such objects. We restrict 
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consideration to area and volume, which is the three-dimensional analogue of perimeter 
and area of two-dimensional models. 

One prominent model is self-avoiding surfaces on the cubic lattice, also studied as 
a model of three-dimensional vesicle collapse. We follow the review in |102] (see also 
the references therein) and consider closed orientable surfaces of genus zero, i.e., surfaces 
homeomorphic to a sphere. Numerical studies indicate that the surface generating function 
displays a square-root 7 = —1/2 as the dominant singularity. 

Consider the fixed surface area ensemble with weights proportional to q n , with n the 
volume of the surface. One expects a deflated phase (branched polymer phase) for small 
values of q and an inflated phase (spherical phase) for large values of q. In the deflated 
phase, the mean volume of a surface should grow proportionally to the area m of the surface, 
in the inflated phase the mean volume should grow like m 3 / 2 with the surface. Numerical 
simulations suggest a phase transition at q = 1 with exponent = 1. This indicates that 
a typical surface resembles a branched polymer, and a concentrated distribution of volume 
is expected. Note that this behaviour differs from that of the two-dimensional model of 
self-avoiding polygons. 

Even relatively simple subclasses of self-avoiding surfaces such as rectangular boxes [73] 
and plane partition vesicles [50], generalising the two-dimensional models of rectangles and 
Ferrers diagrams, display complicated behaviour. Let p m ^ n denote the number of surfaces of 
area m and volume n and consider the generating function P(x, q) = J2mnPm,nX m q n - For 
rectangular box vesicles, we apparently have P(x, 1) ~ A\ log(l —x)\/(l — x) 3 / 2 as x — > 1~, 
some some constant A > 0, see J73J Eq (35)]. In the fixed surface area ensemble, a linear 
polymer phase < q < 1 is separated from a cubic phase q > 1. At q — 1, we have = 2/3, 
such that typical rectangular boxes are expected to attain a cubic shape. We expect a limit 
distribution which is concentrated. For plane partition vesicles, it is conjectured on the 
basis of numerical simulations (5DJ Sec 4.1.1] that P(x, 1) ~ Aexp(a / (x c — x) 1 ^ 3 ) / (x c — x) 1 , 
where 7 ~ 1.7 at x c = 0.8467(3), for non-vanishing constants A and a. It is expected that 
0=1/2. 

As in the previous subsection, three-dimensional models of punctured vesicles may 
be considered. The above arguments hold, if the exponent satisfies < < 1. A 
corresponding result for punctures of unbounded size can be stated if the critical surface 
area generating function is finite. 

3.10 Summary 

In this section, we described methods to extract asymptotic area laws for polygon models 
on the square lattice, and we applied these to various classes of polygons. Some of the laws 
were found to coincide with those of the (absolute) area under a Brownian excursion and a 
Brownian meander. A combinatorial explanation for the latter result has not been given. 
Is there a simple polygon model with the same area limit law as the area under a Brownian 
bridge? The connection to stochastics deserves further investigation. In particular, it would 
be interesting to identify underlying stochastic processes. For an approach to a number of 
different random combinatorial structures starting from a probabilistic viewpoint, see [82] . 
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Area laws of polygon models in the uniform fixed perimeter ensemble q = 1 have been 
understood in some generality, by an analysis of the singular behaviour of g-functional 
equations about the point (x,q) = (x c , 1). Essentially, the type of singularity of the 
half-perimeter generating function determines the limit law. A refined analysis can be 
done, leading to local limit laws and providing convergence rates. Also, limit distributions 
describing corrections to the asymptotic behaviour can be derived. They seem to coincide 
with distributions arising in models of punctured polygons, see [51] . 

For non-uniform ensembles, concentrated distributions are expected, but general re- 
sults, e.g. for g-functional equations, are lacking. These may be obtained by multivariate 
singularity analysis, see also [21] 165]. 

The underlying structure of g-functional equations appears in a number of other combi- 
natorial models, such as models of two-dimensional directed walks, counted by length and 
area between the walk and the x-axis, models of simply generated trees, counted by the 
number of nodes and path length, and models which appear in the average case analysis 
of algorithms, see [51] [37]. Thus, the above methods and results can be applied to such 
models. In statistical physics, this mainly concerns models of (interacting) directed walks, 
see [18] for a review. There is also an approach to the behaviour of such walks from a 
stochastic viewpoint, see e.g. the review [101] . 

There are exactly solvable polygon models, which do not satisfy an algebraic g-difference 
equation, such as three-choice polygons [H], punctured staircase polygons [45J, prudent 
polygon subclasses [96], and possibly diagonally convex polygons. For a rigorous analysis 
of the above models, it may be necessary to understand g-difference equations with more 
general holonomic solutions, as q approaches unity. 

Focussing on self-avoiding polygons, it might be interesting to analyse whether the 
perimeter generating function of rooted self-avoiding polygons might satisfy an implicit 
equation Eq. (I3"3l) . Asymptotic properties of the area can possibly be studied using stochas- 
tic Loewner evolution [60]. Another open question concerns the area limit law for q ^ 1 
or the perimeter limit law for x ^ x c , where Gaussian behaviour is expected. At present, 
even the simpler question of analyticity of the critical curve x c (q) for < q < 1 is open. 

Most results of this section concerned area limit laws of polygon models. Similarly, 
one can ask for perimeter laws in the fixed area ensemble. Results have been given for 
the uniform ensemble. Generally, Gaussian limit laws are expected away from criticality, 
i.e., away from x = x c . Perimeter laws are more difficult to extract from a g-functional 
equation than area laws. We will however see in the following section that, surprisingly, 
under certain conditions, knowledge of the area limit law can be used to infer the perimeter 
limit law at criticality. 

4 Scaling functions 

From a technical perspective, the focus in the previous section was on the singular be- 
haviour of the single- variable factorial moment generating function gk{x) Eq. (TTJ, and on 
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the associated asymptotic behaviour of their coefficients. This yielded the limiting area 
distribution of some polygon models. 

In this section, we discuss the more general problem of the singular behaviour of the two- 
variable perimeter and area generating function of a polygon model. Near the special point 
(x,q) = (x c , 1), the perimeter and area generating function P(x,q) = Ylm>oP'ni{.q)x rn = 
^2n>o a n(x)q n is expected to be approximated by a scaling function, and the correspond- 
ing coefficient functions p m (q) and a n (x) are expected to be approximated by finite size 
scaling functions. As we will see, scaling functions encapsulate information about the limit 
distributions discussed in the previous section, and thus have a probabilistic interpretation. 

We will give a focussed review, guided by exactly solvable examples, since singularity 
analysis of multivariate generating functions is, in contrast to the one-variable case, not 
very well developed, see [81] for a recent overview. Methods of particular interest to poly- 
gon models concern asymptotic expansions about multicritical points, which are discussed 
for special examples in [80, 5J. Conjectures for the behaviour of polygon models about 
multicritical points arise from the physical theory of tricritical scaling [H], see the review 
[6T] , which has been adapted to polygon models [HI [13] . There are few rigorous results 
about scaling behaviour of polygon models, which we will discuss. This will complement 
the exposition in [17]. See also [121 Ch 9] for the related subject of scaling in percolation. 



4.1 Scaling and finite size scaling 

The half-perimeter and area generating function of a polygon model P(x, q) about (x, q) = 
(x c , 1) is expected to be approximated by a scaling function. This is motivated by the 
following heuristic argument. Assume that the factorial area moment generating functions 
gk(x) Eq. ([I]) have, for values x < x c , a local expansion about x = x c of the form 

where 7^ = {k — 6i)/<p and 9i+\ > Q\. Disregarding questions of analyticity, we argue 

In the above calculation, we replaced P(x,q) by its Taylor series about q = 1, and then 
replaced the Taylor coefficients by their expansion about x = x c . The preceding heuristic 
calculation has, for some polygon models and on a formal level, a rigorous counterpart, 
see the previous section. In the above expression, the r.h.s. depends on series J-i{s) = 
S&>o( — -0* fk,is~ lktl of a single variable of combined argument s = (1 — x/x c )/(l — qY . 
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Restricting to the leading term I = 0, this motivates the following definition. For > 
and x c > 0, we define for numbers s_, s + G [— oo, +oo] the domain 

D(s_, 8+ ) = {(x, q) G (0, oo) x (0, 1) : s_ < (1 - x/x c )/(l - qf < s + )}. 

Definition 4 (Scaling function). For numbers p m ,n with generating function P(x,q) = 
J2mnP™.,nX m q n , let Assumption^ be satisfied. Let < x c < 1 be the radius of convergence 
of P(x, 1). Assume that there exist constants s_,s + G [— oo, +oo] satisfying s_ < s + and 
a function T : (s_, s + ) — ■> R, swc/i i/iai g) satisfies, for real constants 9 and <fi > 0, 

FW( I)g )~(l- g ) J f^j (x,q)^(x c ,l) inD(s„,s + ). (34) 

T/ien, the function T(s) is called an (area) scaling function, and 9 and (ft are called critical 
exponents. 

Remarks, i) In analogy to the one- variable case, the above asymptotic equality means 
that there exists a power series convergent for \x\ < X\ and |g| < q±, where 

x\ > x c and q\ > 1, such that the function P( sin9 \x,q) := P(x,q) — P( re9 \x,q) is asymp- 
totically equal to the r.h.s.. 

ii) Due to the region D(s-, s+) where the limit (x, q) — > (x c , 1) is taken, admissible values 
(x, q) satisfy < q < 1 and < x < x (q), where x (q) = x c (l — s_(l — 5)*), if s_ 7^ —00. 
Thus, in this case, the critical curve x c (q) satisfies x c (q) > x (q) as q approaches unity. 
Note that equality need not hold in general. 

iii) The method of dominant balance was originally applied in order to obtain a defining 
equation for a scaling function .F(s) from a given functional equation of a polygon model. 
This assumes the existence of a scaling function, together with additional analyticity prop- 
erties. See [H EQ E7]. 

iv) For particular examples, an analytic scaling function J-'(s) exists, with an asymptotic 
expansion about infinity, and the area amplitude series F(s) agrees with the asymptotic 
series, see below. 

v) There is an alternative definition of a scaling function [21] by demanding 

in a suited domain, for a function H(t) of argument t = (1 — q)/ (l — x/x c ) 1 ^. Such a scaling 
form is also motivated by the above argument. One may then call such a function Ti(t) a 
perimeter scaling function. If T{s) is a scaling function, then a function T~C(t), satisfying 
Eq. (1351) in a suited domain, is given by 

H(t) = t e F{r*). 

If s_ < and s + = 00, the particular scaling form Eq. (|34p implies a certain asymptotic 
behaviour of the critical area generating function and of the half-perimeter generating 
function. The following lemma is a consequence of Definition HI 
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Lemma 3. Let the assumptions of Definition^ be satisfied. 

i) If s + = oo and if the scaling function J-{s) has the asymptotic behaviour 

F{s) ~ f s"y° (s - oo), 
then 7q = — 4, and £/ie half-perimeter generating function P(x, 1) satisfies 

P^\x, 1) ~ /o(l - x/a; c )W (I / x c ). 

7/s_ < and i/ie scaling function J-{s) has the asymptotic behaviour 

F(s) ~ /i s Q ° (s \ 0), 
i/ien «o — 0, and the critical area generating function P(x c , q) satisfies 

p( sin °\x c ,q)~h (l-q) e 



□ 



A sufficient condition for equality of the area amplitude series and the scaling function 
is stated in the following lemma, which is an extension of Lemma [3J 

Lemma 4. Let the assumptions of Definition^ be satisfied. 

i) Assume that the relation Eq. (3$ ) remains valid under arbitrary differentiation w.r.t. 
q. If s + = oo, if the scaling function J-{s) has an asymptotic expansion 

k>0 

and if an according asymptotic expansion is true for arbitrary derivatives, then the 
following statements hold. 

a) The exponent 7*. is, for k 6 N , given by 

k-e 



b ) The scaling function J-{s) determines the asymptotic behaviour of the factorial 
area moment generating functions via 



1 d k 

k\W P{x ' q) 



(sing) 



Jh {x / X c ). 



(1 - x/x c )^ 
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ii) Assume that the relation Eq. (34\ ) remains valid under arbitrary differentiation w.r. 
to x. If s_ < ; and if the scaling function T{s) has an asymptotic expansion 

F{s) ~ ^(-l) fc /^ ( S \ 0), 

fc>0 

and if an according asymptotic expansion is true for arbitrary derivatives, then the 
following statements hold. 

a) The exponent otk is, for k G N , given by a>k = k. 

b ) The scaling function determines the asymptotic behaviour of the factorial perime- 
ter moment generating functions at x = x c via 

( 1 d k \ {smg) x- k hi. 

where /3 k = k(fi — 9 . 

□ 

Remarks. Lemma H] states conditions under which the area amplitude series coincides 
with the scaling function. Given these conditions, the scaling function also determines the 
perimeter law of the polygon model at criticality. 

In the one-variable case, the singular behaviour of a generating function translates, 
under suitable assumptions, to the asymptotic behaviour of its coefficients. We sketch 
the analogous situation for the asymptotic behaviour of a generating function involving a 
scaling function. 

Definition 5 (Finite size scaling function). For numbers p m>n with generating function 
P(x,q) = J2mnPm,nX m q n , let Assumption^ be satisfied. Let < x c < 1 be the radius of 
convergence of the generating function P(x, 1). 

i) Assume that there exist a number t + G (0, oo] and a function f : [0,t + ] — > M., such 
that the perimeter coefficient function asymptotically satisfies, for real constants 70 
and <fi > 0, 

[x m }P(x, q) ~ x^m^fim^il - q)) (q, m) -> (1, 00), 

where the limit is taken for m a positive integer and for real q, such that m l ^{l — q) £ 
[0,t+]. Then, the function f(t) is called a finite size (perimeter) scaling function. 

ii) Assume that there exist constants t_ G [— 00, 0), t + G (0, 00], and a function h : 
[t_ , t + ] —>■ M., such that the area coefficient function asymptotically satisfies, for real 
constants fio an d <fi > 0, 

[q n ]P(x, q) ~ n A_1 /i(n^(l - x/x c )) (x, n) -> (x c , 00), 

where the limit is taken for n a positive integer and real x, such that n^(l — x/x c ) G 
[£_,£+]. Then, the function h(t) is called a finite size (area) scaling function. 



40 



Remarks, i) The following heuristic calculation motivates the expectation that a finite 
size scaling function approximates the coefficient function. For the perimeter coefficient 
function, assume that the exponents of the factorial area moment generating functions 
are of the special form 7^ = (k — Q)/<j). We argue 



00 



w]p{z, q ) - K] g(-i)* _ - q y 



fc=0 

00 „ 



fc=0 

In the above expression, the r.h.s. depends on a function f(t) of a single variable of 
combined argument t = m 1//< ^ ) (l — q). 

For the area coefficient function, we assume that (3k = kef) — 9 and argue as above, 

[q n ]P(x, q) « [<f ] £(-1)*—** (1 - 

^n*-^(-l) fc -^(^(l-x/x c )) fe . 

In the above expression, the r.h.s. depends on a function h(t) of a single variable of 
combined argument t = rv{l — x/x c ). 

ii) The above argument suggests that a scaling function and a finite size scaling function 
may be related by a Laplace transformation. A comparison with Eq. (1201) leads one to 
expect that finite size scaling functions are moment generating functions of the limit laws 
of area and perimeter. 

Hi) Sufficient conditions under which knowledge of a scaling function implies the existence 
of a finite size scaling function have been given for the finite size area scaling function [T3] 
using Darboux's theorem. 

A scaling function describes the leading singular behaviour of the generating function 
P(x,q) in some region about (x, q) = (x c , 1). A particular form of subsequent correction 
terms has been argued for at the beginning of the section. 

Definition 6 (Correction-to-scaling functions). For numbers p m ,n with generating function 
P(x,q) = J2mnPm,nX m q n , let Assumption^ be satisfied. Let < x c < 1 be the radius of 
convergence of the generating function P(x, 1). Assume that there exist constants s_, s + G 
[—00, +00] satisfying s_ < s + , and functions T\ : (s_, s + ) — > R for I 6 No, such that the 
generating function P(x, q) satisfies, for real constants <f> > and 61, where 9i + i > 9\, 

P^(x,q) ~ Y^-zfTi (iT^v) (Xc ' 1} m D ^ s ^- 

Then, the function Tq^s) is a scaling function, and for I < 1, the functions Ti(s) are called 
correction-to-scaling functions. 
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Remarks, i) In the above context, the symbol ~ denotes a (generalised) asymptotic 
expansion (see also [HOI Ch 1]): Let (Gk(x)) keNo be a sequence of (multivariate) functions 
satisfying for all k the estimate Gk+i{x) = o(Gk{x)) as x — > a? c in some prescribed region. 
For a function G(a?), we then write G(x) ~ Xlfclo^^) as x ~~ > ^ f° r a ^ ra we have 
Cf^) = Y2=l G k{x) + 0(G n (cc)) asx ^ x c . 

ii) The previous section yielded effective methods for obtaining area amplitude functions. 
These are candidates for correction-to-scaling functions, see also [87J. 



4.2 Squares and rectangles 

We consider the models of squares and rectangles, whose scaling behaviour can be explicitly 
computed. Their half-perimeter and area generating function can be written as a single 
sum, to which the Euler-MacLaurin summation formula [80, Ch 8] can be applied. We 
first discuss squares. 

Fact 5 (cf [4"9l Thm 2.4]). ForO < x,q < 1, the generating function P{x, q) = Ylm=o xm( l m2 ^ 4 
of squares, counted by half-perimeter and area, is given by 

P(x,q) = ^=^=F[ '^IL \+l + R(x,q), 




y/\logq\ \ a/I log? I / 2 
with J-"(s) = ^/vr e s2 erfc(s), where the remainder term R(x,q) is bounded by 

\R(x,q)\ < -| logx|. 
o 



□ 



Remarks, i) The remainder term differs from that in [521 Thm 2.4], where it was estimated 
by an integral with lower bound one instead of zero [IHJ Eq. (46)]. 

ii) With x c = 1, S- = and s + = oo, the function J-"(s) is a scaling function according to the 
above definition. The remainder term is uniformly bounded in any rectangle [xo, 1) x [q , 1) 
for < xq, qo < 1, and so the approximation is uniform in this rectangle. 
Hi) The generating function P(x,q) satisfies the quadratic g-difference equation P(x,q) = 
l-\-xq 1 ^P(q l l 2 x, q). Using the methods of the previous section, the area amplitude series of 
the model can be derived. It coincides with the above scaling function ^{s) . This particular 
form is expected, since the distribution of area is concentrated, p(x) = 5(x — 1/4), compare 
also with Ferrers diagrams. 

iv) It has not been studied whether the scaling region can be extended to values x > 1 
near (x, q) = (1,1). It can be checked that the scaling function J-'(s) also determines 
the asymptotic behaviour of the perimeter moment generating functions, via its expansion 
about the origin. As expected, they indicate a concentrated distribution. 

The half-perimeter and area generating function of rectangles is given by 



x(qx) 7 



q' x 

=1 s=l r=l ^ 
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We have P(x, 1) = x 2 / (1 — x) 2 , and it can be shown that P(l, q) ~ _l£l(Lji as g y \ 



9 



sec 



[851119]. The latter result implies that a scaling form as in Definition HI with s_ < 0, does 
not exist for rectangles. We have the following result. 

Fact 6 ([HJ Thm 3.4]). For < q < 1 and < qx < 1, the generating function P(x, q) of 
rectangles satisfies 

P(x, q) = tt— — r f 7T~~~T - LerchPhi (qx, 1, ) J + R(x, q), 

I log g | V|logx| V I log ^ I / / 

with the Lerch Phi-function LerchPhi (z, a, v) = Y^=o ( v + n ) a > w ^ ere the remainder term 
R(x, q) is bounded by 

x 2 q (I |logx|^ x 2 q | log g 

\R{x,q)\ < z - + 



1 — qx \2 6 / (1 — qx) 2 6 

□ 

Remarks, i) The theorem implies that, for every go £ (0, 1), the function (1 — qx) 2 P(x, q) 
is uniformly approximated for points (x, q) satisfying q < q < 1 and < x < x c (q), where 
x c (q) = 1/q is the critical curve. 

ii) Rectangles cannot have a scaling function ^(s) as in Definition H] with s_ < 0, since 
the area generating function diverges with a logarithmic singularity. This is reflected in 
the above approximation. 

Hi) It has not been studied whether the area moments or the perimeter moments at criti- 
cality can be extracted from the above approximation. 

iv) The relation of the above approximation to the area amplitude series of rectangles of 
the previous section, F(s) = Ei(s 2 ) e s , is not understood. Interestingly, the expansion of 
F(s) about s = resembles a logarithmic divergence. It is not clear whether its expansion 
at the origin is related to the asymptotic behaviour of the perimeter moment generating 
functions. 



4.3 Ferrers diagrams 

The singularity diagram of Ferrers diagrams is special, since the value x c (l) := lim^i x c (q) 
does not coincide with the radius of convergence x c of the half-perimeter generating function 
P(x, 1). (The function q t— > x c (q) is continuous on (0,1], as may be inferred from the 
exact solution.) Thus, there are two special points in the singularity diagram, namely 
(x, q) = (x c , 1) and (x, q) = (x c (l), 1). Scaling behaviour about the latter point has 
apparently not been studied, see also [85] . 

About the former point (x, q) = (x c , 1), scaling behaviour is expected. The area ampli- 
tude series F(s) of Ferrers diagrams is given by the entire function 

F(s) = y|erfc (V2s) e 2 * 2 . 
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A numerical analysis indicates that its Taylor coefficients about s = coincide with the 
perimeter moment amplitudes at criticality, which characterise a concentrated distribu- 
tion. There is no singularity of F(s) on the negative real axis at any finite value of s, in 
accordance with the fact that the critical line at q = 1 extends above x = x c . 

It is not known whether a scaling function exists for Ferrers diagrams, or whether it 
would coincide with the amplitude generating function, see also the recent discussion [501 
Sec 2.3]. An rigorous study may be possible, by first rewriting the half-perimeter and 
area generating function as a contour integral. A further analysis then reveals a saddle 
point coalescing with the integration boundary at criticality. For such phenomena, uniform 
asymptotic expansions can be obtained by Bleistein's method [HQl Ch 9.9]. The approach 
proposed above is similar to that for the staircase model [83J in the following subsection. 



4.4 Staircase polygons 

For staircase polygons, counted by width, height, and area with associated variables x, y, q, 
the existence of an area scaling function has been proved. The derivation starts from an 
exact expression for the generating function, which has then been written as a complex 
contour integral. About the point (x,q) = (x c , 1), this led to a saddle-point evaluation 
with the effect of two coalescing saddles. 

Fact 7 (cf [831 Thm 5.3]). Consider < x, y, q < 1 such that the generating function 
P(x, y, q) of staircase polygons, counted by width, height and area, is convergent. Set 
q = e~ e for e > 0. Then, as e\0, we have 

P{x,y,q) = (- | — -+ 



Ai(ae" 2 / 3 ) 

uniformly in x,y, where a = a(x,y) satisfies the implicit equation 

\c?l* = log(x) |° g( ^~^ + 2U 2 (z m - V~d) -2U 2 (z m + Vd), 
o log(z m + Vd) 

where z m = (1 + y — x)/2 and d = z^ — y, and Li 2 (t) = — f* log ^~"' > du is the Eul 
dilogarithm. □ 

Remarks, i) The characterisation of a 3 / 2 given in [831 Eq (4.21)] has been used. 

ii) The above approximation defines an area scaling function. For x = y and x c = 1/4, we 

obtain the approximation [831 Eq (1.14)] 

P, x ,)^l + 4-^ Al ' (44/3(1/4 -^ 2/3) 



er 
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as (x,q) — > (x c , 1) within the region of convergence of P(x,q). It follows by comparison 
that the area amplitude series coincides with the area scaling function. 
Hi) An area amplitude series for the anisotropic model has been given in [56], by a suitable 
refinement of the method of dominant balance. 

iv) It is expected that the perimeter law at x = x c may be inferred from the Taylor ex- 
pansion of the scaling function T{s) at s — 0. A closed form for the moment generating 
function or the probability density has not been given. The right tail of the distribution 
has been analysed via the asymptotic behaviour of the moments [5?] [55] . See also the next 
subsection. 

v) The above expression gives the singular behaviour of P(x,q) as q approaches unity, 
uniformly in x, y. Restricting to x = y, it describes the singular behaviour along the line 
q = 1 for < x < x c . In the compact percolation picture, this line describes compact per- 
colation below criticality. Perimeter limit laws away from criticality may be inferred along 
the above lines. (Asymptotic expansions which are uniform in an additional parameter 
appear also for solutions of differential equations near singular points [80] .) 

vi) By analytic continuation, it follows that the critical curve x c (q) for P(x,x,q) coincides 
near q = 1 with the upper boundary curve x$(q) = (1 — s_(l — g) 2 / 3 )/4 of the scaling 
domain, where the value s_ is determined by the singularity of smallest modulus of the 
scaling function on the negative real axis, hence by the first zero of the Airy function. This 
leads to a simple pole singularity in the generating function, which describes the branched 
polymer phase close to q = 1. 



4.5 Self-avoiding polygons 

In the previous section, a conjecture for the limit distribution of area for self-avoiding 
polygons and rooted slef-avoiding polygons was stated. We further explain the underlying 
numerical analysis, following [91]. [9?. 93J. The numerically established form Eq. (132]) implies 
for the area moment generating functions for k ^ 1 singular behaviour of the form 

with critical point x c = 0.14368062927(2) and 7& = 3k/2 — 3/2, where the numbers fk are 
related to the amplitudes A k in Eq. fl32]) by 

A i. h 



For k — 1, we have 7^ = 0, and a logarithmic singularity is expected, gi(x) ~ /ilog(l — 
x/x c ), with fx = A\. Similar to Conjectured] this leads to a corresponding conjecture 
for the area amplitude series of self-avoiding polygons. If the area amplitude series was 
a scaling function, we would expect that it also describes the limit law of perimeter at 
criticality x = x c , via its expansion about the origin. (Interestingly, these moments are 
related to the moments of the Airy distribution of negative order, see [93] 134"].) This 
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prediction was confirmed in [93], up to numerical accuracy, for the first ten perimeter 
moments. Also, the crossover behaviour to the branched polymer phase has been found 
to be consistent with the corresponding scaling function prediction. As was argued in the 
previous subsection, the critical curve x c (q) close to unity should coincide with the upper 
boundary curve xo(q) = x c (l — s_(l — q) 2 ^ 3 ), where the point s_ is related to the first 
zero of the Airy function on the negative real axis, s_ = —0.2608637(5). The latter two 
observations support the following conjecture. 

Conjecture 3 ([871 [93]). Let p m ,n denote the number of self-avoiding polygons of half- 
perimeter m and area n, with generating function P(x,q) = 'Yli mn Pm,nX m q n ■ Let x c = 
0.14368062927(2) be the radius of convergence of the half-perimeter generating function 
P(x, 1). Assume that 

^p m ,n ~ A x- m m- 5/2 (m oo), 

n 

2 

where A is estimated by A = 0.09940174(4). Let the number s_ be such that (4A ) 3 7rs_ 
coincides with the zero of the Airy function on the negative real axis of smallest modulus. 
We have s_ = -0.2608637(5). 

i) For rooted self-avoiding polygons with half-perimeter and area generating function 
p( r \x,q) = x4^P(x, q) , the conjectured form of a scaling function JF^s) : (s_,oo) — ► 
K as in Definition ^ is 

^ (r)(s) = 27d7 logAi ( (4Ao)§7rs )' 

with critical exponents 6 = 1/3 and = 2/3. 

ii) The conjectured form of a scaling function J-(s) : (s_,oo) — > K for self-avoiding 
polygons is obtained by integration, 

Ha) = log Ai (( 4 ^o) f vrs) + ^(1 - q) log(l - q), (36) 
with critical exponents 6 = 1 and <fi = 2/3. 

Remarks, i) The above conjecture is essentially based on the conjecture of the previous 
section that both staircase polygons and rooted self-avoiding polygons have, up to normal- 
isation constants, the same limiting distribution of area in the uniform ensemble q = 1. 
For a numerical investigation of the implications of the scaling function conjecture, see the 
preceding discussion. 

ii) A field-theoretical justification of the above conjecture has been proposed [IE]. Also, 
the values of A\ = l/(87r) and the prefactor l/(127r) in Eq. (1361) have been predicted using 
field-theoretic methods [15], see also the discussion in [93] . 
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4.6 Models in higher dimensions 

Only very few models of vesicles have been studied in three dimensions. For the simple 
model of cubes, the scaling behaviour in the perimeter-area ensemble is the same as for 
squares |4*9l Thm 2.4]. The scaling form in the area- volume ensemble has been given [4"9l 
Thm 2.8]. The asymptotic behaviour of rectangular box vesicles has been studied to some 
extent [73]. Explicit expressions for scaling functions have not been derived. 

4.7 Open questions 

The mathematical problem of this section concerns the local behaviour of multivariate 
generating functions about non-isolated singularities. If such behaviour is known, it may, 
under appropriate conditions, be used to infer asymptotic properties such as limit distri- 
butions. Along lines of the same singular behaviour in the singularity diagram, expressions 
uniform in the parameters are expected. This may lead to Gaussian limit laws [37]. Parts 
of the theory of such asymptotic expansions have been developed using methods of sev- 
eral complex variables [81]. The case of several coalescing lines of different singularities 
is more difficult. Non-Gaussian limit laws are expected, and this case is subject to recent 
mathematical research [8~T] . 

Our approach is motivated by certain models of statistical physics. It relies on the 
observation that the singular behaviour of their generating function is described by a 
scaling function. There are major open questions concerning scaling functions. On a 
conceptual level, the transfer problem [35] should be studied in more detail, i.e., conditions 
under which the existence of a scaling function implies the existence of the finite-size scaling 
function. Also, conditions have to be derived such that limit laws can be extracted from 
scaling functions. This is related to the question when can an asymptotic relation be 
differentiated. Real analytic methods, in conjunction with monotonicity properties of the 
generating function, might prove useful [80] . 

For particular examples, such as models satisfying a linear g-difference equation or di- 
rected convex polygons, scaling functions may be extracted explicitly. It would be interest- 
ing to prove scaling behaviour for classes of polygon models from their defining functional 
equation. Furthermore, the staircase polygon result indicates that some generating func- 
tions may have in fact asymptotic expansions for q /* 1, which are valid uniformly in the 
perimeter variable (i.e., not only in the limit x /* x c ). Such expansions would yield scaling 
functions and correction-to-scaling functions, thereby extending the formal results of the 
previous section. This might be worked out for specific models, at least in the relevant 
example of staircase polygons. 
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